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PREFACE 


The conclusions given in this report are based upon the "Engineering 
Manual (PM100-1)" method for culaculation of "protection factors". Since 
an error analysis is not presently available, the conclusions should be 
regarded as tentative, pending the development of such an analysis. In 
addition, a redistribution of fallout end/or changes in the 7-ray spectrum 
emitted by the fallout may introduce further uncertainties into these con¬ 
clusions. 
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ABSTRACT 
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This is Volume II of four separately bound volumes that report the research 
completed in fulfillment of Office of Civil Defense Work Unit No. 3233B, "Radio¬ 
logical Recovery Requirements, Structures, and Operations Research." This describes 
six supporting studies all previously reported to the Office of Civil Defense in 
research memoranda. Volume I describes the general aspects of the investigations 
and presents the conclusions and recommendations. The abstract for each of the 
volumes is presented on the following pages. 
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ABSTRACT FOR VOLUME I 


This study examines the application of decontamination strategies to 
extensive urban areas. Urban areas of various sizes (from a few acres to an 
interconnected system involving hundreds of acres) are examined with regard to 
decontaminating vital sections and their connecting links. Feasible creation of 
decontaminated "islands" or marshalling areas is determined. The nature and 
scope of command and control system elements required for conducting effective 
decontamination in practical situations are defined together with the preattack 
and postattack data required by such a system. Several models of detector 
response to gamma radiation developed during the course of the project are briefly 
discussed. 
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ABSTRACT FOR VOLUME II 


Volume II contains six supporting studies. Described are a number of models 
for determining the cost and effectiveness of decontaminating .municipal areas. 

One study examines the system components of a command and control system for 
municipal decontamination. These studies cover the following subjects: 

1. A Feasibility Study of the Application of Analog Computers to the Analysis 
of Decontamination : This study describes the research completed on the applica¬ 
tion of analog computers to the analysis of decontamination. A simplified analog 
model of the effect of a single decontamination effort on a single plane of 
contamination is explained and described. Sample analog records of output from a 
prototype of this model are presented. The design of a more elaborate model 
requiring substantially more analog equipment is described, and applications are 
indicated. 

2. A Circular Model for Approximating Gamma Ray Intensity at a Single 
Detector Location : An approximate procedure for determining the gamma ray 
intensity at a point, due to fallout radiation, is investigated. A reasonably 
accurate model using circular annulus geometry was developed which includes the 
effects of gamma ray attenuation, build-up, backscatter, and skyshine. This 
model was based on NBS Monograph U2. 

A single annulus model is described in order to show clearly the basic 
premises upon which the full "circular" model is based. The development of the 
multiple annulus model is treated in a manner to show how each of the parameters 
involved is handled. An example demonstrating the use of the model is included. 

3. A Square-Grid Model for Approximating Gamma Ray Intensity at a Sl..gle 
Detector Location : This study describes a simple and practical procedure for 
determining the approximate fallout gamma ray intensity at a point as a function 
of the geometry of the contributing planes and related shielding. The procedure 
described employs a "square-grid" technique for modeling the contaminated planes. 
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This model is shown to be reasonably accurate and easily adapted to practical 
situations with the help of scaled overlays for city maps to be used in conjunc¬ 
tion with the tables and graphs developed in the research. An example analysis 
is carried out to illustrate the method. 

4. A Point Source Model for Approximating Gamma Ray Intensity at a Single 
Detector location (The Equivalent Planes Method ): A method for analyzing the 
dose contributions from plane areas of fallout is developed; each contaminated 
plane is treated as a weighted point source of radiation intensity. Only the 
area, the location of the center, the eccentricity of the approximating rectangle, 
and the intervening shielding need be considered. Limits are developed for the 
eccentricity of an area, centered over the detector, which can be represented by 
a square of the same irea with an error of less than ten percent in dose rate 
contribution. Limits are also developed for the area of off-center contaminated 
areas such that the product of area tim<j>s the dose rate contribution per unit 
area is within ten percent of the true contribution from the area. 

A sample analysis is given. 

A step by step procedure based on this model, called "The Equivalent Planes 
Method," is presented. A small booklet which could be used as an "Equivalent 
Planes Method" workbook (together with some work sheets) is included in an envelope 
at the back of this volume. 

5. A FORTRAN Program for Decontamination Analysis : This study describes a 
debugged and tested FORTRAN computer program to compute the effectiveness parameters 
used to analyze municipal decontamination. The program was written in FORTRAN 

64 tc he used on large scale computers such as the CDC 3600, 

6 . The Nature and Scope of Command and Control System Elements Required lor 
Conducting Effective Decontamination in Municipalities } This study serves to 
determine the nature and scope of command and control system elements which are 
required to effect practical municipal decontamination. The preattack and 
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postattack data requirements f or decontamination are specified, ard the essential 


components of 
The influence 


an information system for decontamination are identified and related, 
of direct weapons effects on the decontamination system is examined. 




ABSTRACT FOR VOLUME III 


Volume III contains the cost and effectiveness data for decontamination 
analyses of sixteen sites and facilities in San Jose, California. Costs are 
measured in team-hours of effort. Decontamination effectiveness is measured In 
terms of fractions of dose-rate remaining at specified detector locations and 
fractions of dose remaining for persons who perform functions requiring specified 
daily activity patterns at the sites and facilities chosen. 


x 




ABSTRACT FOR VOLUME IV 


Volume IV contains the cost and effectiveness data for decontamination 
analyses of twelve sites and facilities in Detroit, Michigan. Costs are measured 
in team-hours of effort. Decontamination effectiveness is measured in terms of 
fractions of dose-rate remaining at specified detector locations and fractions of 
dose remaining for persons who perform functions requiring specified daily activity 
patterns at the sites and facilities chosen. 
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Radiological Recovery Requirements, Structures, and Operations Re s earch 
Volume ~I, Development of Analytical. Computer, and Systems Models it. Support 
of Decontamination Analysis 

I. INTRODUCTION 

This volume reports six studies performed under Office of Civil Defense 
Subtask 3233B, Radiological Recovery Requirements, Structures, and 
Operations Research. It is addressed to technical personnel concerned with the 
planning of postattack recovery operations as summarized and described in Volume 
I of this report. All of the studies presented in Volume II are concerned with 
the development of tools which can be used to examine the effectiveness and 
costs of decontamination when applied to accelerating recovery of an activity in 
a postattack environment. 

When this work began, it seemed that large complexes comprised of several 
buildings could not be analyzed quickly or efficiently without modifying the 
analysis techniques that were previously developed for single-facility shielding 
analysis. Thus, to meet what then appeared to be a requirement for efficient 
decontamination analysis, a number of models were developed for approximating 
gamma ray intensity at a point due to complex contaminated plane configurations. 
Each of these models is discussed in a separate appendix (Appendices A through 
D) in this volume. 

Paralleling the development of the analytical models and the analog computer 
model, two computer programs written in FORTRAN for the CDC 3690 were completed 
and debugged to perform most of the computation required to analyze decontamination 
operations as applied to several sites and facilities selected from San Jose and 
Detroit. The first of these programs (Reference 1) was developed and debugged 
under another contract. This program computes the plane-by-plane contributions to 
intensity at a specified detector location. The second of these programs computes 
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the effectiveness parameters of the individual decontamination analyses. This 
program is described in Appendix E. 

Lastly, the nature and scope of the command and control system elements 
required to affect practical municipal decontamination are defined. The pre¬ 
attack and postattack data requirements are identified and related to the system 
as a whole. The influence of direct weapons effects on the decontamination system 
are examined. This work is described in Appendix F. The six supporting studies 
included in this volume are: 

1 . Appendix A: A Feasibility Study of the Application of Analog Computers 
to the Analysis of Decontamination 

2. Appendix B: A Circular Model for Approximating Gamma Ray Intensity at 
a Single Detector Location 

3. Appendix C: A Square-Grid Model for Approximating Gamma Ray Intensity 
at a Single Eetector Location 

4. Appendix D: A Point-Source Model and the Equivalent Planes Method for 
Approximating Gamma Ray Intensity at a Single Detector Location 

5. Appendix E: A FORTRAN Program for Decontamination Analysis 

6 . Appendix F: The Nature and Scope of Command and Control System Elements 
Required for Conducting Effective Decontamination in Municipalities 
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A Feasibility Study of the Application of Analog 
Computers to the Analysis of Decontamination 


Note; The material in this Appendix was originally submitted to USNRDL 
as Research Memorandum RM-OU-214-2*. 


* J. T. Ryan. A Feasibility Study of the Application of Analog Computers 
to the Analysis of Decontamination . RM-OU-214-2. Durham, North 
Carolina: Research Triangle Institute, Operations Research and 
Economics Division, 1 April 1965 
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Appendix A 


A Feasibility Study of the Application of 
Analog Computers to the Analysis of Decontamination 

I. INTRODUCTION 

To the author's knowledge, only digital computer programs are presently used 
to calculate the protection factor associated with specific detector locations. 
These programs, with some modifications, can be used to calculate the individual 
intensity contributions from each of the contributing contaminated planes. Thus, 
these programs can be effectively used to estimate the value of decontaminating 
one or more of these planes of contamination. These programs are very detailed. 
The inputs are an itemized description of all of the planes of contamination and 
of the intervening shielding. 

The over-all time required to prepare such inputs is usually very long, and 
thus reduces the number of facilities or activities which can be considered in 
the analysis of decontamination for a large municipal area. This deficiency led 
to the investigation of an analog model that uses continuous variables for the 
environmental parameters as well as tho parameters associated with the decontami¬ 
nation operations. Analog models are characteristically smaller in scope than 
digital models, less minute in detail, have a shorter running time, and can be 
more easily changed. 

This study is a result of a very brief effort to determine the applicability 
of analog computers to the analysis of decontamination. This appendix presents 
a detailed, non-technical description of the analog computer study. The initial 
pages are devoted to a description and discussion of analog computer?. The sub¬ 
sequent pages describe both the limited decontamination simulation fitted to the 
available hardware at the Research Triangle Institute and a more elaborate model 
which requires substantially more equipment. Sample runs are included. The final 
section of the appendix presents conclusions and recommendations for future work. 
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Section VI 
ware components 


io a technical description of the mathematical equations and hard- 
(both existing and proposed) of t-he decontamination models. 



II. DISCUSSION OF ANALOG COMPUTERS* 


A. General 

The history of computers is, in fact, two histories. Host authors (References 
A-l and A-2 for example) recognize that there are very basic factors which dis¬ 
tinguish between the family of digital computers and the family of analog com¬ 
puters. Usually the abacus and the slide rule are used to exemplify a simple 
computer from each of these families respectively. For the purpose of this report, 
however, it is not necessary to elaborate on the philosophy or historical divergence 
of the two families of computers. It is only necessary to discuss the one primary 
factor which distinguishes between analog and digital coraputers--the way in which 
data is handled in each machine. 

In the digital computer, data is comprised of discrete numbers (represented 
by digits, magnetic senses, bit positions—or in the case of the abacus—bead 
positions, etc.) whereas in the analog computer data is represented as a continuous 
variable (represented by wheel positions, electric voltages—or in the case of the 
elide rule-slide bar positions, etc.). 

The general-purpose electronic analog computer is the type of computer 
considered in this report. In order to evaluate the features of this type of 
analog computer, it is helpful also to consider electronic digital computers. 

The following several paragraphs extracted from Reference A-3 very clearly present 
the pertinent differences between analog and digital models as they affect the 
analysis and conclusions presented in this report: 


*Readers already familiar witn analog computers should omit this section. Readers 
who would like to learn more about the use of analog computers are referred tc 
Reference A-l. 

For a detailed (technical) comparison of digital and analog computers, see 
Reference A-2, Chapter 2. 




"Digital computers are composed of a large number of multiple inter¬ 
connected bistable devices. These on-off components along with trans¬ 
ferring devices form the basis of the digital computer capability-- 
they enable formation and storage of numbers, simple arithmetic opera¬ 
tions, and two-valued logic operations. Any function which can be 
expressed or approximated as a sequence of these operations can be per¬ 
formed on the digital computer. 

"Analog computers are not as simply composed but consist of many 
components with different functions or operating modes. Their capa¬ 
bility depends on their having a large and varied collection of com¬ 
ponents to draw on. Electronic circuits provide this required variety, 
including summers, integrators, multipliers, function generators, 
switches, comparators, and many simple circuits which may be built 
quickly from resistors and capacitors to provide special functions. 

"Mechanizi.-.g an analog computer model consists in connecting these 
various components together in the same way as the equations are formed 
cj as the physical system is envisioned. The model is built in parallel; 
for each operation in the system there is a component in the model that 
perf^; ms the same operation. If there are ten multiplications in the 
system, then there must be ten multipliers in the model, which will all 
be operating at once. This is just the opposite of digital computers, 
which are organized ip series. That is, every operation in the system 
is performed by one component. If there are ten multiplications in the 
system, then the digital model will perform one multiplication at a 
time in the programmed sequence. 

"One consequence of parallel organization is that increases in 
system size or complexity bring a corresponding increase in the number 
of components required in the model. With series organization, an 
increase in system size or complexity results in a longer running time 
for the model. 

"A second consequence of the parallel organization is that if the 
model is interrupted at any time during the run, all variables will be 
found to be at the value corresponding to the one the system would 
have at that time. That is to say, there is an isomorphism between 
time in the system and time in the model. This is not the case with 
series organization as in digital computer models, where the model can 
be interrupted only at certain specific times to be meaningful. 

"Computing time on an analog computer model is usually very short. 
While the basic switching time of the digital computer may be one 
thousand to ten thousand times faster than the response tine of an 
analog computer component, the number of operations required in the 
digital computer for any complicated function is so large that, in 
fact, the total operation may not be faster than the analog computer 
component. Thus, speed in performing the total operation coupled 
with parallel organization allows the time variable on the analog to 
be scaled as a fraction or multiple of real time. Running times on 
the analog are usually scaled to range between 10 seconds and 2 
minutes. For example, an analog model of a valve which opens in 2 
seconds might be scaled to have 1 second of real time equal 10 seconds 
of computer time, thereby slowing do”n the operation for each observa¬ 
tion. A battle model in which the time of flight to target is 5 hours 
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might be scaled to have 1 hour of real time equal 1 second on the 
computer, vastly speeding up the time required for a run. 

"Analog computer models require a considerably more contracted 
treatment of large-scale problems but provide a more tractable model. 

That is, large volumes of data, readily handled by digital computers, 
are impossible on the analog unless they can be treated in an aggre¬ 
gated manner; but changes can be made easily, while operating, by 
simply resetting a dial. 

'There is considerable difference between the form of the output 
on the two computers. The usual output of a digital computer is a 
tabulation of the results. Analog output is usually in a graphic 
form, whether as a chart showing the variation of a number of para¬ 
meters with time or as a plot of one variable against another. This 
graphical output is adapted very well to problems where a physical 
understanding of the interactions is desired, since the effect of 
variations can be immediately seen. 

"In summary then, an analog computer model with its treatment of 
data and variables quite aggregated would be expected tc be considerably 
smaller in scope than a digital computer model. Operational changes 
could be simply made and input data easily varied. Output would be 
graphical and running time short, so that results of input changes 
would be immediately available and easily interpreted. This is in 
contrast with present digital models which are detailed, longer run¬ 
ning, less readily changed, and considerably broader in scope." 

B. Analog Computers and Decontamination Analyses 

Some comments augmenting and tying the above remarks to the problem of 
analyzing decontamination operations are appropriate here. 

The "parallel organization" feature described above is particularly useful in 
analyzing the effects of decontaminating a number of independent contaminated 
surfaces. Instead of having to compute the effects of the separate operation 
serially, one can build the model so as to show, as a function of (scaled) real 
time, the effect of the sum of these decontamination operations, where they could 
be performed simultaneously, or sequenced in any way whatsoever. 

Furthermore, different differential equations can be used to govern the 
efficiency of the individual decontamination operations where the parameters 
describing these differential equations can be varied manually at run time. 

This, of couise, would be impossible to do with a digital computer. The ability 
to easily solve very complex differential equations (not easily soluable by paper 




and pencil methods) is one of the chief capabilities of analog computers. 


The "graphical output" feature can, of course, be performed on a digital 
computer with the appropriate peripheral equipment, but the ability to respond at 
computer running time to the information being graphed is impossible with normal 
digital equipment. Thus, analog computers can be a useful tool for training 
personnel responsible for large-scale municipal decontamination. This will be 
discussed in more detail later in this appendix. 




III. DESCRIPTION OF I-IUTED DECONTAMINATION SIMULATION 


A. Introduction 

The limited decontamination analog model was developed to explore the possi¬ 
bilities arid determine the application of analog computers to decontamination 
analysis. As such, this model does not provide realistic answers to decontamina¬ 
tion problems and only roughly simulates any real situation. Very gross treatment 
of the details was used ir. order to simulate the main features of a single decon¬ 
tamination operation on a small sized computer. 

B. The Hardware Used 

The limited analog model described in this secLion was fitted to a basic 
installation of the model 3400 Desk Top Donr.er Analog Computer (Figure A-l). The 
basic model 3400 computer contains ten amplifiers, usable as summers or integrators. 
The Research Triangle Institute installation also includes a Donner Model 3430 
Problem Board and a Donner Model 3073 Potentiometer Strip. The problem board 
(shown in Figure A-2) is used to interconnect components for the solution of 
particular problems and the potentiometers are used to adjust constant and 
parametric coefficients in the equations which describe particular problems. The 
output from the analog model is graphed using an Offner Dynograph Amplifier- 
Recorder Model 542. The recorder is also shown on che right in Figure A-l. A 
more detailed description of the actual components used n ti.is model is con¬ 
tained in section VI. 

C. Th e Limited Decontamination Mode l 

In the limited decontamination analog, each contaminated plane which con¬ 
tributes to the intensity at a detector location is approximated by a point source 
at the centroid of the contaminated plane. 

Figure A-3 is a schematic showing the relationships between a single detector 
and four contaminated planes. 
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Figure A-2 

A Close-up of the Problem Board 
Set up to Analyze a Simple 
Decontamination Operation 
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Figure A-3 


Schematic Diagram of Simplified Decontamination Operations 
with Four Contaminated Planes. 




The contributions, C., to radiation intensity from ground level sources of 


contamination are computed using Equation A-l. 



(A-l) 


where * the area in square feet of the i th contaminated plane, 

r^ = the diPtance in feet from the centroid of the i^ contaminated plane 
to the de'.ector, and 

= the shielding factor associated with the i** 1 contaminated plane. 

Her 2 represents the fraction of the total intensity (received at the detector) 
whi :h comes from the i th plane. The shielding factor is a dimensionless number 
which attenuates the contribution according to the shielding between the detector 
and the i th plane. (Note that 0 < < 1). 

The fraction of the total intensity received at the detector which comes from 


the roof of the building in which the detector is centrally located is found by 
*- 

Equation A-2 


2 

c = n In (1 + ^ ) s 
h 


(A-2) 


where r « the radius of a circular roof, with an area equal to the area of the 
actual roof over the detector, 
h = the height of the roof above the detector, and 
= the shielding factor associated with the roof. 

Equation A-2 can also be used to determine the contribution to the intensity at an 
unshielded detector (out-of-doors) from the plane directly below the detector. 

Here = 1 and h is a height corresponding to the ground roughness factor. 


*See page 743 of Reference A-7 for a derivation of Equation A-2. 

See Reference A-5 for an explanation of how the ground roughness factor affects 
the "effective" height of the detector. 
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In the simplified model both -j S and n In (1 + -j ) S are treated as single 
r 1 h i 

constant inputs. Any variations in these individual parameters can be simulated 
by appropriately changing the capacitors which represent these aggregate variables 
before running the model. 

Decontamination is governed (for each individual plane) by Equation A-3. 

* * -K.E 

F i = F i + (1 " V 6 (A " 3) 

th 

where = the fraction of fallout remaining on the i contaminated plane after 
it has been decontaminated. 

F* = the limiting fraction of fallout remaining on the i** 1 contaminated 
plane after infinite decontamination. 

K. = a constant associated with a given method and the physical nature of 
the i fc k contaminated plane. 

E. - the effort (usually measured in man-hours) which is applied to the 
i th contaminated plane. 

Here, both 1 - F* and E. are treated as single variable inputs. Both of 

these variables can be varied by changing settings of the potentiometer which 
it? assigned to the particular variable parameter. This can be done while the 
model is being run. 

Table A-I is a summary of the equations used in the limited model of decon¬ 
tamination. 

While the model is being run, the individual contributions are s”mmed and 
the radiation intensity" at the detector is output as a function of time. The 
total dose than an individual would receive at the detector location is also out¬ 
put as a function of time. The dose is measured over the time interval from 
to t 2 - This is simply given by 

**See Reference A-6 for a derivation of Equation A-3. 
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operations, H(t) will be a monotonically decreasing step function. 
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TABLE A-I 


System Equations £fc Simplified Model 

1. Intensity contribution from the i h ground level contaminated plane: 

A i 

C. - 3. where 

1 'i 1 

A^ * the areas in square feet of i t ^ 1 contaminated plane. 

r. - distance in feet from centroid of fallout on i tn contaminated 
plane to the detector location. 

S^ = shielding factor associated i^ contaminated plane. 

2. Intensity contribution from roof: 

2 

C r = H In (1 + ) S^ where 

h 

r * the radius of a circular roof of equal area, 
h « height of roof (above detector). 

S = shielding factor associated with roof. 


3. Decontamination efficiency: 


F i = F* + (i - F*) e" K i E i where 

F^^ = fraction of fallout remaining after decontaminating i^ plane 

Fj^ = fraction of fallout which cannot be removed. 

= constant associated with i 1 ”* 1 surface and the method used to 
decontaminate it. 

E^ = the amount of effort applied to decontaminate the i^ contami 
nated plane. (Usually measured in man-hours of effort). 


For the simplified model the following parts of the above equations were set to 
single constant parameters; 


i A i 

k i 4 s i 


2. K r B H In (1 + ~s-) S 
2 h 2 r 


The following parts of the above equation were set to single variable parameters 

. ■> r,* - 


A-13 




IV. SAMPLE RUNS WITH PROTOTYPE MODEL 

Several test runs were made. Seven f these runs are included here to 
illustrate the nature and format of the answers obtainable. The situation 

simulated in all seven runs was the decontamination of a single contaminated 

-1 2 

plane. For simplicity, no fallout build-up function was used. The t ’ decay 
curve was approximated as the sum of two exponential decay curves (i.e., 
t =* e K 1 C + e K 2 fc ). A fuller discussion of this approximation is contained 
in Chapter VI. 

In all of these sairple runs, no attempt was made to quantitatively scale 
dose-rate and total dose. The time variable was also left arbitrary. Thus, two 
inches of horizontal distance on the graph outputs might correspond to two weeks, 
two hours, etc. The actual graph outputs from these seven runs are presented 
as Figures A-4 through A-10. Time runs from right to left on these graphs. For 
Figures A-4 through A-9, decontamination is assumed to take n small finite time. 

In Figure A-10 decontamination is instantaneous. 

Figure A-4 shows the effect of decontaminating a single plan of contamination 
where 90% of the fallout material has been removed. Figures A-5 and A-6 show 
the effect of decontaminating the same plans of contamination where the decontami¬ 
nation efforts expended are only 75% and 50% of the effort expended in the opera¬ 
tion simulated by the run which produced Figure 4. (Still assuming the percent 
of removable fallout equals 90%). 

Figure A-7 shows the effect of decontaminating the same plane where only 
50% of the fallout can be rimoved and a "full" effort (the same effort is n«sumed 
for Figure A-4) is expended. Figure A-8 shows the effect of decontaminating 
this plane where only 50% of "full" effort is expended. 

For all of the Figures A-4 through A-8, both total dose ..nd dose-rate sre 
shown where total do-*e begi is to accumulate at time t^ and stops accurw lat ing at 
time t^- Note that the vertical scale for total dose in Figures A-7 and A-8 is 
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exactly half of the scale used for total dose in Figures A-4, A-5, and A-6. 

In Figure A-9, total dose begins to accumulate at where no build-up 
function is assumed. The total dose continues to accumulate throughout the graph 
on Figure A-9. Here, the horizontal scale for total dofj was again cut to half of 
thrt used in Figure A-8. The same fraction of removable fallout and effort 
expended are used in Figures A-8, A-9, and A-10. 

Figure A-10 shows how the graphs would look if decontamination took place 
instantaneously. Otherwise Figure A-10 is the same as Figure A-9. 
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Figure A-4 

Analog Record of Dose Rate and Total Dose as a Function of Time 
for Case I 



instantaneous time of arrival of fallout (no buildup), 
time when decontamination begins. 

t^) * interval of time for which total dose is calculated. 


Thus, 



t < t 
t 1 < t < C. 


t > t 


2 
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Figure A-5 

Analog Record of Dose Rate and Total Dose as a Function of Ti 
for Case II 


Legend 

t = time of arrival of fallout (no buildup). 
a 

t * time when decontamination begins, 
d 

(t^, t 2 ) = interval of time for which total dose is calculated, 


Thus, 














Figure A-6 

Analog Record of Dose Rate and Total Dose as a Function of Time 
for Case III 


Legend 

t = time of arrival of fallout (no buildup). 
a 

tj = time when decontamination begins. 

^1’ ^2^ = ^ n ^ erv<i ^ f° r which total dose is calculated. 



V. 


t < t 

1 


t < t 


2 






Figure A-7 

Analog Record of Dose Rate and Total Dose as a Function of Time 
for Case IV 

id 

t a “ time of arrival of fallout (no buildup). 
t^ = time when decontamination begins. 

(t,, t )= interval of time for which total dose is calculated. 








Figure A-8 

Analog Record of Dose Rate and Total Dose as a Function of Time 
for Case V 


Legend 

t “ time of arrival of fallout (no buildup). 


(t^, t^) = interval of time for which total dose is calculated. 
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Figure A-9 

Analog Record of Dose Rate and Total Dose as a Function of Time 
for Case VI 


Legend 

t & = time of arrival of fallout (no buildup). 

t = time when decontamination begins, 
d 

Total dose is calculated from t to any time t . 


Thus, 


Dose = 
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V. AN ELABORATE MODEL 


Although quite adequate to demonstrate feasibility, the limited analog 
model described in Section III is not suitable for complex decontamination prob¬ 
lems. A more elaborate model with less aggregation of variables and a capacity 
for handling a larger number of contaminated planes is necessary. The basic 
system equations, however, remain the same for the larger analog model described 
in this section. 

As in the model already implemented and described in Section III, decontami¬ 
nation is measured by a "graph" of dose-rate (radiation intensity) at some detector 
location of interest. Total dose is measured for selected time intervals at the 
same point. Complex buildup functions can be used to simulate the initial stages 
of fallout arrival. The number of contaminated planes which can be considered for 
a given detector, as well as the number of detector locations, is in theory, 
unlimited. A practical limit on the numbers of analog components would limit the 
number of contaminated planes to 20 and the number of detector locations to three. 
By appropriately using relays and switching circuits, the total dose received by 
a person changing environments instantaneously (such as moving from indoorc to 

outdoors) could be simulated. 

-1 2 

All functions (log, t ' , etc.) can be quite accurately simulated by using 
suitably adapted function generators. This will improve the precision of the 
simulation. 

Furthermore, a quick-response capability could be included to serve as an 
analysis and training tool. This is accomplished by using manual over-rides 
which permit the operator to apply decontamination efforts on the basis ol graphs 
which portray previous decontamination efforts during the same run. 

A technical description and flow diagram of the more elaborate model is 
included in Chapter VI. 
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VI. A TECHNICAL DESCRIPTION OF THE ANALOG DECONTAMINATION MODELS 


A. Calculations to be Performed 

The calculations to be performed by the computer are expressed in Equations 
A-5 and A-6: 



D(t) = / R(t)dt = f v(t; t 1 ,t 2 ) R(t)dt (A-6) 


where the dymbols are defined as follows: 

R(t) = dose rate in roentgens per hour, 

D(tj * dose (normalized units) in roentgens, 

n = number of contaminated planes considered, 

A^ = area of i^ contaminated plane in square feet, 

r^ = distance from detector to centroid of i t * 1 contaminated plane in feet, 
S i = shielding factor for i tl1 location (C < < 1), 

K, - constant for i C ^ contaminated plane and a particular decontamination 

1 method, ^ 

E^ = effort applied to i^ contaminated plane in man hours/1000 sq ft, 

Cg = effort coefficient * 1 - e ^i^i, 

i 

C = fraction of removable mass for i*"* 1 contaminated plane, 
m. 

x i = decontamination time constant, 

p = effective radius of surrounding contaminated plane (roof or detector 
surface), 

h “ effective height of surrounding contaminated plane (roof or detector 
surface), 

■ shielding factor for surrounding contaminated plane (roof or detector 
" s surface), 
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i 


“ integration limits for dose D(t). 

The function u(t-t- ) is a unit step function which has the values: 
i 


u(t-t ) 
L i 


t < t < « 
i “ 


The function vCt^it,,^) has the values: 

v(t > t 1 ,t 2 ) =0 0 < t < t 1 

=1 ti < t < t 2 

= 2 t ? < t < =o 

A functional diagram showing the calcjlations to be performed is shown in Figure 
A-ll. Only one of the n identical channels (each channel corresponding to a 
single contaminated plane) is shown. For n areas, the calculations represented 
within the dotted outline will be repeated n times. 


B. Ran ge of Constants, Parameters, and Variables 

The units of the various quantities and the expected ranges of these quantitic 


are listed below: 



SYMBOL 

UNITS 

RANGE 

R(t) 

r/hr 

- 

t 

hrs 

IU -336 (K 00) 

n 

dimensionless 

1-10 

A 

ft 2 

100-100,u00 

r 

ft 

10-500 

S 

dimensionless 

0-1 


dimensionless 

0-1 

K 

2 

ft ,'raan-hour 

.1-5 
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SYMBOL 

UNITS 

RANGE 

G 

man hrs/1000 ft 2 

.01-1 

C 

dimensionless 

0-1 

m 

T 

hr" 1 

10-0.1 

\ 

hr 

10-1000 

P 

ft 

10-200 

h 

ft 

1-100 

S 

dimensionless 

0-1 

s 

*2 

hr 

10-1000 

fc l 

hr 

10-1000 

D(t) 

roentgens 

- 


For computational purposes the dose rate R(t) and dose P(t) are measured in 
normalized units. To convert to roentgens per hour or roentgens, one point on the 
dose rate curve muse be hand calculated, and the proper scale factor determined. 

C. Computer Scaling 

1. Choice of Time Scale 

The time range of major interest is from 24 to 336 hours. For convenience, 
the computer t' ^ scale, t c> may be chosen so that one hour of real time, 
t„, correspond* to .5 second of computer time. Time zero on the computer 
will correspond to time 24 hours in real time. Thus: 

t c = ~ (t r - 24) seconds, (A-6) 

where t £ s: computer tim- in seconds, and 

t r = real time in hours. 

2. Dose Rate Scale 

At t r = 24 hrr, the dose rate R(t) will be assir t ed the valus 100 in 
arbitrary units. The computer scale factor for R(t) will be taken as 1 
volt/unit. To convert from arbitrary units to roentgens per hour, the scale 
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factor is determined by hand calculation of the dose rate at t^ » 24 hrs. 

The dose D(t) will be scaled at .01 volt per arbitrary unit to 
prevent the integrator from overloading with long integration times. 

3. Coefficient Scaling 

For adjusting values of the various parameters that remain constant 
during a particular computer run, the computer potentiometers will be used. 

In order to simplify the equipment requirements, it is convenient to group 
certain parameters into dimensionless ratios where possible. These parameter 
groups will t! be set into the computer as potentiometer settings. Addi¬ 
tional coefficients are defined as follows: 

A i 

C = —- range .01 - 10 

r i r- 1 

1 + fiy I range 0-20 

h J 

For potentiometers with dial divisions from 0 to 100, the following scale 
factors are assigned to the parameters that are to be set on potentiometer 
dials: 

5c* i e Factor 

(units per dial division) 

.1 

.2 

.01 

.01 

.01 
.01 


Coefficient 


V 


Range 

.01-10 

0-20 

0-1 


0-1 

0-1 
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It should be noted that with additional equipment complexity, each value sucn 
as A^, etc. may be entered separately as a dial setting, and the necessary 
multiplications and divisions accomplished by the computer. 

4. Dynamic Range Requirements 

In some cases, it would be desirable to simulate a time period longer 
than 15 days. In increasing the time period of a single computer run, diffi¬ 
culties may be encountered because of tue small signal levels that will be 

- 1.2 

present. For example, the function t decreases by a factor of approxi¬ 
mately 100 between 24 hours and 1100 hours. If the computer output were 
initially set at 100 volts, the output at 1100 hours (550 seconds computer 
time) would be approximately 1 volt if no decontamination had taken place. 

This level is not too small; however, the signal level from each of the i 
channels must be considerably less than 1 volt (with 10 channels giving equal 
contributions, the signal level would be .1 volt). On certain computers, 
the noise level prevents accurate operation at these low voltage levels. 

In order to avoid the inaccuracies inherent in using low signal levels, 
the best approach would appear to be to break the problem into time periods 
of interest. If the period of interest were from 10 to 20 days, the computer 
time would be scaled such that the problem starts at t^ = 240 hrs instead of 
24 as discussed previously. 

D. Computer Diagram 

A computer diagram to perform the necessary calculations is shown in Figure 

A-13 where the symbols are defined in Figure A-12. A function generator generates 
-1 2 

a voltage proportional to t . This voltage is adjusted to a value analogous to 

the contribution of each of the contaminated planes by the potentiometers C , S , 

r i 1 

C fl , and S g . The voltages analogous to the radiation received from each area are 
summed in a summing amplifier, and the resultant voltage is recorded on a strip 
chart recorder. 
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Figure A-12 

Table of Flow Chart Symbols and Operations 
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To simulate decontamination of a given area, the relay Ry is closed. A 
time lag of the form e t/T i is introduced by operational amplifier #7, and the 
original voltage is reduced by the proper amount by subtraction in amplifier #1. 

The amount of voltage subtracted depends upon the effort coefficient and the per¬ 
cent of removable material as set on the potentiometers C_ and C . 

^i m i 

To calculate the dose from the dose rate, an integrator is energized by relay 
Ry B. This relay remains closed for a time period t ^ as set on the sequencing 

function generator. 

Note that the computer diagram shows only one of n identical setups that 
must be patched into the computer. 

E. Function Generators 

-1 2 

To simulate a function varying with time as t , several methods are avail¬ 
able. A conventional diode function generator with a linear time function input 
may be used. Another method is to use an x-y plotter as a function generator, 

supplying a linear time function as a drive for the x axis, and reading out a 
-1 2 

voltage proportional to t * on the y axis by means of a curve following probe. 

-1 2 

The simplest technique, however, appears to be that of approximating the t 
function by a sum of exponentially decaying functions. For example, with three 
exponentials, the t curve may be fitted by the sura at six points on the curve. 
The exponential functions are generated very easily on the analog computer by 

setting up an initial charge on a capacitor, and discharging the capacitor slowly 

-1 2 

through a resistor. Figure A-14 shows the approximation of the t ' function 
in the range 24-350 hrs by the sum of two exponentials as follows: 

t' 1 ' 2 * ^ [.86 e- 03U(t - 24) + .14 .-«»»(«-»>] 

The computer diagram for accomplishing this approximation would appear as shown 
in Figure A-15 (.5 sec computer time B 1 hr real time). 
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Figure A-15 
Approximation of t 


Function 


The functions u (t-t^) and v (t, t^, are generated by a sequencing device 
which may be a set of motor driven cams. Switches are closed by the cams at the 
desired time, and the proper relays are closed as shown in Figure A-12. Provision 


must be made for convenient adjustment of the times t^, 
driven cams, n+1 cams will be required. 



contaminated areas, the following analog equip- 


F. Equipment Requirements 

To set up a situation with 10 
ment will be required: 

Quantity 

34 

52 
( 1 ) 
or (3) 


lt< T> 

Operational Amplifiers 

Dial Setting Potentiometers 

* 1.2 

Function Generator t 

Operational Amplifiers (for 
-1 2 

approximation of t ' function) 
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Quantity 

11 

1 

1 


Item 

Relays 

Sequence Function Generator 
(11 Channel) 

Dual Channel Recorder 


A 40 amplifier analog computer should be sufficient for the 10 area simulation. 

The sequence function generator may be constructed In the laboratory if a suitable 
type cannot be purchased. To buy this equipment would cost about $40,000. Total 
rental costs (including set-up costs) would be under $1,000 for a one year period. 
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VII. CONCLUSIONS AND RECOMMENDATIONS 


A. Conclusions 

The basic conclusion of this study is that the analog computer provides a 
potentially useful and quite feasible tool for simulating and analyzing decontami» 
nation operations . 

The development and implementation of this first analog model of decontamina¬ 
tion has pinpointed the areas of difficulty and the direction in which further work 

appears most promising. The most difficult problem encountered was accurately 

- 1.2 

approximating the t decay rate. This problem is considered in detail in 
Chapter VI; since a high-accuracy analytic function generator would be very costly 

(about $1000), some form of approximation involving the sum of exponential functions 
-K t 

of the form K^e 2 is a practical compromise solution. 

Even though a number of the problems associated with the simple mod^l using 
only one plane of contamination were a result of complex scaling necessary to keep 
the range of variables within machine limits, future models, to be useful, must 
include even more detail and so must use even more equipment. Thus, scaling will 
become an increasingly difficult problem. 

B. Future Work 

Any future work must be undertaken in the light of the major difficulty found 
in using the limited a-.i&log model--the large amount of equipment required for an 
adequate simulation. Nonetheless it is recommended that two major directions be 
explored: 

!• The development of a larger more detailed simulation, but still with 
simplified aggregation of variables where possible, end 
2. The development of a quick response capability to serve as a ready 
analysis and training tool. 
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GLOSSARY OF TERMS 


A = Intercept value used in the analytical approximation of 

L (X,o>) functions. 

d « Vertical height of the detector above the source plane. 


D q » Total detector response at a standard reference position. 

= Gamma dose rate received at a detector from the annular 
region. 


D 


kj 


Gamma dose rate received at a detector from the annular 
sector of the annulus. 


D/D * Dose reduction factor, 

o 


h 


kj 


Height of contaminated annulus above the ground or reference 
plane. 


L (X) * Total detector response for infinite plane isotropic source. 

L (X,cu) * Generalized expression for the geometry factor. 

L c (X,to) ■ Geometry factor describing detector response due to circular 

plane sources of fallout radiation. 

M “ Fractional area of the k** 1 annular region covered by the 

J contaminated j th annular sector. 

* Total number of sectors in the k tn annulus. 

S - Slope value used in the analytical approximation of L (X,o) 

functions. 

0 ) ■ Solid angle fraction. 

X * Effeccive mass thickness. 

* Height of k** 1 annular region above the detector. 



A Circular Model for Approximating Gamma Ray 
Intensity at a Detector Location 

I. INTRODUCTION 

The purpose of this appendix is to show the initial results of recent studies 
directed toward the development of rather simple and practical procedures for 
determining the gamma ray intensity at a point due to nuclear weapon induced 
fallout. The consensus of those involved in this research effort was that if 
such an approach proved feasible, it would result in a considerable reduction 
in both labor and computer costs for those cases where apprc /mate results are 
considered adequate. In addition, the development of such an analytical tool 
would be of considerable value in simplifying those tasks associated with the 
ensuing analysis of decontamination effectiveness in municipal areas. This 
research is a first step in the development of several alternative approaches to 
a model for field analysis of fallout decontamination effectiveness. 

In an effort to assure a reasonable degree of accuracy in the analytical model 
presented in the following discussions, such pertinent factors as gamma ray attenu¬ 
ation. build-up, backscatter, and skyshine, are incorporated. Another basic 
premise closely adhered to in connection with the development of this model is 
that it be built around the familiar results contained in NBS Monograph 42 
(Reference B-l). 

The general notational scheme presented here, together with selected elementary 
geometric considerations and the information contained in Reference B-l, provides 
the basis for what is called a "circular model". This model is utilized to furnish 
values of the so-called "reduction factor," D/D o , due to various planes of con¬ 
tamination at different heights with or without the inclusion of barrier shielding. 








In particular, these planes of contamination are divided into annular sectors wit., 
the detector located In each instance at the origin of the polar coordinate system. 

The approach taken in describing this circular model is that of developing 
initially the detector response for the single contaminated annulus case (Section 
II) following with increasingly more complex annular geometries (Section III) and 
terminating with a brief discussion on barrier sheilding considerations (Section IV). 

Each of the annular geometric configurations utilized in the analytical 
development of the model required the extensive use of the L (X) and the L(X,w) 
functions developed in Reference B-l. It was found during the course of the 
investigation that the L (X,tu) functions could be approximated over their range 
of greatest utility by simple analytical expressions (Section VI). In turn, such 
expressions can be quickly and easily evaluated with the aid of a desk calculator 
and standard mathematical tables, or efficiently stored in a computer memory without 
the necessity of lengthy subroutines. 

The procedure for treating attenuation through barriers makes effective use 
of the L (X) function and the (X,o)) geometry factor, through judicious choices 
and combinations. It should be noted that considerable refinement of these 
procedures and associated calculations may be necessary for situations near the 
detector location at the origin of the coordinate system. 
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II. SINGLE ANNULUS MODEL 


The single annulus is chosen initially since it will permit the development 
of the basic premises from which more complex geometries can later be described. 


A. Complete Un iformly Contaminated Annulus 

In order to enpr'-sb the detector response for the case of a single uniformly 
contaminated annulus (see Figure B-l), it is necessary to first define the solid 
angle tractions and subtended at the detector by the contaminated 

annulus. 


1. Solid Angle Expressions 


From Figure B-l we define first the solid angle fraction 
subtended by angle Os^ as 

cd, * 1 - cos , (See Section V for derivation of solid 

a a angle expressions' 


where 



Similarly, 

*(k«). • 1 • C0S a (k+l)a 


(k-t-l)a 


K(k+l) 


(B-l) 


(B-2) 


(B-3) 


(B-4) 
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2. Detector Response 


The detector response for the single annulus can be computed by 
appropriate differencing of the geometry factor L c (X,u)) (see Figure B-2 
and Section VI) for the values of ^ and This differencing can be 

illustrated in terms of the detector response fov the uniformly contaminated 
annulus in Figure B-l as follows: 

£ - MX) [l c - 1*. (X,^)] , (B-5) 

where D q is the detector response at a standard reference position*, and 
is the response due to the k tl * annulus. The function L (X) is graphed 
in Reference B-l. 

It should be noted that all distances d from the detector to the 
contaminated plane of interest may be expressed in terms of effective mass 
thickness (X) for the medium between the source and detector. For the case 
of air (at 20°C., 76 cm Hg) as this medium, the following relationship may be 
written (see p. 19, Reference B-l): 

d/X - 13.3 Ft. of air/psf. (B-6) 


The interchangeability of d and X 
analytical expressions is permissable 
indicated in Equation B-6. 

B. Sinale Annular Sector Case 

The similar sector shown in Figure B-l 


as the penetration variable in these 
with the use of che conversion factor 

L represent* a fraction 


,!ia—i 

i 2 


(8*7) 


a 

B is defined here as the detector response at an unprotected position in an infinite 
homogeneous air medium located three feet above e hypothetical infinite plance source. 
Thu source it considered to be tf the seme character et the radiation fallout 
uniformly distributed on the contamination pian«(s) being investigated. 
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Figure B-2 


Ceo»etry Factor (L^ <X,w)) Depicting Detector Response 































of the total annulus. Consequently, the detector response due to the annular 
sector In the annulus is simply 


Ski - 

D 


Mj 1 (X) 


^(k+Ua^c 


(B-8) 


For the case of N annular sectors comprising all the contaminated areas 
of the k** 1 annulus, the detector response is given as follows: 


d £2 

j-i 


L (X) 


L c < X ’>+l)a> ' L c (X ’“ka> 
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uniformly contaminated annular sector 




III. MULTIPLE ANNULAR REGION MODEL 


The following discussion outlines the expansion of the previously developed 
concepts foi a single uniformly contaminated annulus to that of the more frequently 
encountered multiple annular region case. 

The notation utilized in developing mathematical expressions for this model 
will relate to the physical description shown in Figure B-3. The use of the 
subscript, j , which refers to the (here) nonexistent annular sectors, is not 
required, but is retained for consistency in notation. 



First it is necessary to determine the detector response due to the 
k th annular region. This quantity is given by the expression 



* k is the height of the k annular region above the detector. 



Figure B-3 

Multipit Annular hagions couflgurat -n Shoving 
Unifornly Cootaainated Regions of Arbitrary Haight 
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2. Contaminated Regions of Arbitrary Height 


With the k*"* 1 annular region detector response defined, we next explore 
a condition which has considerable practical usage when attempting to deter¬ 
mine detector response in urban areas where radioactive fallout can be 
considered uniformly distributed. It is desired to compute the detector 
response to N annular uniformly contaminated regions of arbitrary height. 

The detector response from a set of such annular regions is given by 



Here it is assumed that each annular region is not effectively "shadow" 
shielded by another. A more generalized analysis concerning mutual shielding 
and barrier attenuation will be given in Section IV. 

As a means of illustrating the summing of the detector responses 
associated with equivalent and uniformly contaminated annular planes of 
different h^ values, the following example is included. 

From the information provided in Figure B-4, an expression for the total 
detector response can be written as follows: 


ID 
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Figure B-4 

Detector Reeponae Due to Uniformly Contaminated 
Fiance with Different Values 
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k=3 
L D 
D_ _ k= l. 


D 0l < h l + 3 > 


D. (d-h.) D, (d-h„) D (d-h ) 

= L ^ r + v^ ) + V^ 


(B-12) 


where D ■ D * D = D 
°1 °2 °3 ° 


B. Equivalent And.Uniformly Contaminate d Annular Sectors for a Set of Annular 
Regions■ 

The detector response which is attributed to the annular sector of the 
k th annular region may be defined by 


- M L (X) 
a 


- \ 

L \ W(k+l) Z & + (d'h kj )y 

- L 1 "‘V ^ 

^Cka ) 2 + (d-h kj ) 2 / 


(B-13) 


where is again the fractional area of the k annular region covered by the 
contaminated j th annular sector. 

It is then apparent that the detector response “ of the k tk annulus is 
found by summing the expressions as follows where we let ? k be the total 

number of sectors in the k C ^ annulus: 

J‘ p v n 


w ° k r i 

]T, if : 2^ V w [ 1 ' (d 'V“<kH).> ■ L c “’V “k.)j 


(B-14) 
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By knowing the detector response for the k annulus, it is now possible to 
express the detector response for any number of such annular regionj by 
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IV. BARRIER ATTENUATION 


One of the most common situations encountered in the analysis of contaminated 
regions located in core areas of large municipalities is the frequent occurrence 
of physical barriers between the detector and a contaminated plane of interest. 

As an example of a rather simplified approach to the barrier problem, we will 
consider here the case of two annular sectors of separate annular regions. The 
heights of the respective annular sectors are assumed to be such that one is 
completely shielded by the other. This particular case is illustrated in Figures 
B-5 and B-6. 

The gamma radiation dose at the detector from a shielded contaminated 
annular sector is essentially due to: (1) direct gamma radiation passing through 
the barrier; and (2) scattered gamma radiation transported by various interaction 
processes in air to the detector location. Here for simplicity the actual radiation 
dose from a shielded contaminated annular sector will be approximated by 
(see Figures B-6 and B-7) 

- M j L [ “c < x - “(U + l).> - L c <*• “k.> ] ' 

where X is the effective distance through the barrier. This approximation will 
not be found to properly treat the scattered radiation contribution to the detector. 
The functions L (X), L. (d,u ), and L £ (X,cd) limit the detector response to a 
circular area of "ntform gamma radiation intensity whose center is on the perpen¬ 
dicular from the detector to the pi ne and separated by either a thickness of air d 
or barrier thickness X (see figure B-7). A numerical example of these considera¬ 
tions is given in Section VII. 


B-15 



Figure B-5 


Two Annular Sectors of Separate Annular Regions - 
One Sector Shielded From the Detector 



Figure B-6 

Detector Response From a Shielded Sector 
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Figure B-7 

Detector Geometry Associated With (X,u) and 
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V. DERIVATION OF SOLID ANGLE FRACTION EXPRESSION 

The purpose of introducing the concept of solid angle in this memorandum is 
due primarily to the fact that the chosen direction for viewing the geometrical 
configurations associated with the circular model is not limited to that of a 
plane. As a result the more generalized notion of angle, i.e., solid angle, is 
utilized in a somewhat modified manner. 

One of the most common means of envisioning a solid angle is to draw a cone 
frow a point as the apex. Then the ratio of the area intercepted by that cone on 
any sphere centered at the apex to the square of the radius of the sphere is a 
measure of the solid angle 0 subtended by the cone. When measured in this 
manner the unit of solid angle is called the steradian. Since the total area of 
2 

a spherical surface is 4jtR , the solid angle subtended by any surface which 
completely surrounds a point cannot exceed 4 jt steradians. It is then also 
apparent that the maximum solid angle subtended on one side of a plane surface 
is 2* steradians. On the basis of this premise, the unit of angular measurement 
oo adopted in this memorandum wa3 chosen to be 2n steradians. This unit is 
consistent with that found in Reference B-l, where ci is designated as the "solid 
angle fraction" in an effort to avoid confusion in angular terminology. 

By using the parameters indicated in Figure B-8, the derivation of the basic 
solid angle fraction expression found in Section II-A-1 of this appendix can be 
obtained 

With the geometry defined as illustrated above, an expression for the area 
intercepted on the spherical surface can be written as 
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Figure B-8 

Solid Angle Ge one try 
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2jt (R sin a) R da 


2 

-2jtR co3 a 


2«R (1 - cos a^) 


From Che basic definition of a solid angle, the following relationship 


“k. ■ 7 ■ 2 ” <l • oos “k.* 


Since our unit of solid angle is 2jt steradians, Equation B-18 must 
be rewritten in the following manner: 


“ka' ®k« 


where i3 the so-called "solid angle fraction.” 


B-21 



VI. APPROXIMATE ANALYTICAL FIT FOR L (X,o>) FUNCTION 


It was discovered during the course of developing the circular model that an 
analytical fit over a rather extensive range of values for the L c (X,to) curve 
found on page 41 of Reference B-l is possible by employing an expression of the 
general form: 


D/D o - L (X,o>) - (B S e M 

The parameter S is obtained from the previously referenced curves by 
applying the following familiar analytical relationship. 

In L (X,uO - In L (X,u>. ) 

5 a .r, . „ -■ — ——L 

In - In 


(B-20) 


(B-21) 


The parameter A is then calculated by solving Equation B-20 for A as 
follows: 

a « IblL i2UiJLJL-§JLflL .a 
A x 


(B-22) 


The best analytical fit for the L c (X.co) curve was achieved by obtain¬ 
ing the L fi (X.uj) values given on page 41 of Reference B-l corresponding to 
and (^2 values of 0.05 and 0.2 respectively and solving for the intercept 
at co ■ 0.2 for each X value selected. The numerical results of these 
calculations are presented in Tabic B-I, and the resultant curve is available 
as Figure B-2 in this appendix. 

It is apparent when comparing the L (X,ai) curves developed here with those 

appearing in Reference B-l that the accuracy of the L (X,oj) values obtained from 

Equation B-20 deteriorate rapidly due to the lack of convergence as o>->0. This 

is particularly true for those values of co > 0.3 and u> > 0.5 associated with the 

L (X, ) function. It should also bi? noted that there are no restrictions associated 
c 

with this model that would prohibit the exclusive use of exact values L £ (X, ) as 

B-2 3 











VII. EXAMPLE 


In an effort to demonstrate the type of situations best handled by the 
circular model, and at the same time effectively summarize the analytical 
methods presented in this appendix, the following example is included. 


The problem is that of determining the c.etector response (reduction factor) 
due to the contaminated annular sectors (A, It, C) in Figure B-9. For the 
purpose of illustration, an value of 0.2 has been arbitrarily selected for 
each of the three planes of contamination to be investigated. 


Contaminated Sector A; 

First we define the solid angle fraction subtended by angle by 


= 1 - cos 


(B- 23) 


where 


ia k = 


_1QQ_ 


r (k) + (d) 


|f (50) + (100)' 


(B-24) 


Therefore, = 0.107. 

We can now convert the distance in air from the source plane into an 
effective mass thickness value by employing the conversion factor provided 
in Equation B-6 as follows: 

X • .tr ' /p ' s ' f * 7 ' 53 ^ ' < B - 25 > 

The detector response for this case is given by the expression 

D a /D q - Mj L (X) [l c (X, a^)l (B-26) 
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By substituting the values of M., X, and cu^ into the previous expression 


B a /D o - 0.2 L (7.53) 


|l c (7.53, 0.107)J 
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Obtaining the (X,uj) value from Figure B-2, the detector response is 
found to be 


D a /D q = 0.2 (0.375) (0.054) * 0.00405 


(B-28) 


Contaminated Sector 


.th 


The expression for the detector response from the j annular sector 
of the contaminated annulus is derived in Section III and appears below 
in only a slightly modified form. 


V D o * M j L < x) 


L ( d-h, , 1 - d-h, \ 

C 1 ) 

' + (d ‘ h k )2 ' 

■ \ ( d -V 1 - _ \ 

V ^(k) 2 + (d-h k ) 2 / 


In this case the effective mass thickness for air is 


13.3 Ft. of 
air psf 


rx 2.26 psf 


(B-29) 


(B-30) 


With the X value and the information provided in Figure B-9, we can 
proceed to find the detector response in the following manner: 


B- 27 



[ L, I 2.26, 1 - 30 _\ 

\ V (60) 2 +"(30)7 

- L c (2.26, 1 - & - - \ ' 

V ft 50?+ (30) 2 y 


0.2 L (2.26) L (2.26, 0.554) - L (2.26, 


0.487), (B-31) 


Here the L (X) and (X,oj) values are obtained from page 33 and page 41 respectively 
of Reference B-l, since the u> values in this case approach the limits of accuracy 
for the approximation method described in Chapter VI. Finally, 


D b /D q - 0.2 (0.57) (0.045) * 0.00513 . (B-32) 


Conta m ina t ed Se ctor C j. 

In finding the detector response for Sector C, a barrier is encountered 
between the detector and the contaminated plane being investigated. From the 
discussion given in Section IV, we can express the detector response for this 
barrier case as follows: 

W V !X> [ L c (X ’“k2> - l c (X '“kl>] ■ (B - 33 > 


The values for uu^ 

and 0^2 are 

found by 



u kl" 1 ’ 

"* °ki' ’ 

' _ 

m 

_ - 0.142, 

(B-34) 




V(60) 2 + 

(100) 2 


u ic2 * 1 

c °* \2 ■ 

1 * 

J2L. 

_ _ - 0.220 . 

(8-35) 




W(*oi ♦ 

(100) 2 
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The X value chosen for the barrier is 75.3, which is equivalent to 0.523 feet 


of concrete. By substitution of these values into the D c /D q expression, we 
obtain 

D c /D q - 0.2 L (75.3) ^ (75.3, 0.22) - L c (75.3, 0.142)j . (B-36) 

By using the values for (X,oj) in Figure B-2, the detector response for Plane 

C is: 

D c /D 0 “ (0.2) (0.038) [ 0.453- 0.285 j - 0.00128 . (B-37) 

Despite certain limitations inherent in the circular model presented here, 
the methods described can produce simple and practical approximate solutions to 
many types of problems requiring knowledge of gamma ray intensity, due to 
weapon-induced fallout, at a single detector location. 
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VIII. SUMMARY AND CONCLUSIONS 


From the information provided in the previous sections, it is apparent that 
the circular model is somewhat restricted in application. The principal limitation 
is that circular or annular planes of contamination are a poor approximation to 
most real situations. However, it may be possible in many cases to reduce non- 
circular geometries of interest to circular regions of equal area. Of course, 
such an innovation must be approached with caution since for certain configurations 
the error introduced in the detector response due to radical changes in the 
original geometry may be quite significant. (Of course, the "Engineering Manual" 
azimuthal sector approach makes these approximations.) It is noted, however, that 
this appendix represents but one of a series of different approaches to the over¬ 
all problem of quick and efficient approximations to the gamma ray intensities due 
to fallout radiation. Other approaches will not have the particular limitations 
noted here. 
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A So • are-Grid Model for Approximating Gamma 

Rav Intensity at a Single Detector location 


I. INTRODUCTION 

A. Purpose 

This appendix describes Che development of a simple and practical procedure 
for estimating fallout gamma ray Intensity at a detector as a function of the 
geometry of the contributing planes of contamination and related shielding. The 
need for such a procedure arises from requirements for analyses of large municipal 
areas with many contaminated planes and detector locations scattered throughout the 
area. Conceptually, the "circular-model" procedure described in Appendix B is 
simple to use, but accuracy suffers when the contaminated planes and barriers do 
not have approximately polar symmetry about the detector location. A square-grid 
model can accurately approximate real situations more easily sir.^e most planes of 
contamination (and most buildings) are rectangular. Further, the square-grid 
model can be applied easily to an arbitrary detector location with the use of 
scaled map overlays. 

B. Some Analysis Considerations 

A problem encountered during the development of the square-grid model was 
the increased difficulty in computing the solid angles subtended at the detector 
by the square areas Adjunct to this problem is the subsequent determination of 
the radiation at a given detector location using Spencer's L-functions (Reference 
C-l). A simple sector-type weighting as was used in the circular model is not 
possible with the square-grid model. 

The solid angles subtended at a given detector location by a square area can 
be computed rigorously. For example, the solid angle o> subtended by any square 
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can be computed (Reference C-l) by the equation 


“ " * £V< € J. ,7i l> “ < >2 T(< 2”>2 ) 4 V (l 3 , ’ > 3 ) " V ( VV ] 


(C-l) 


where. 


<i * 
< x i ,y i ) t " 


*i “ 
T ‘i “ 


T <vV 


+1 if corner i is in quadrant 1 or 3 
i-l if corner i is in quadrant 2 or 4 ’ 

coordinates of corner i, 

|y t l / l*J . 

1*1 / \\\ 

perpendicular distance from the detector to the plane 
containing the plane of contamination under analysis, 
and directly above the origin, and 

solid angle of a rectangle with height z and corners at 

(x i .y i ), (“vV* <- x i»-yi>’ <V y i)- 


A full discussion of Equation C-l and the variables involved is given in Reference 
C-l. Unfortunately, when computing the solid angles for the situations considered 
in this appendix, it was found that the t( e, tj) curves of Reference C-l were not 
sufficient for computing the solid angle to. This is because differences must be 
taken between numbers of approximately the same magnitude. The r functions 
thus were computed from the analytic expression 


'<«•>>) 
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C. The Model 

The general program which wa* followed for the development of a model tor 
approximating fallout gamma ray intcnsitiea at a detector is now described. The 
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model is primarily aimed at the approximate estimation, in the field, of decontamina¬ 
tion effectiveness. High speed computer programs are not practical in these 
circumstances. Of course, building and area characteristics must be known 
approximately. A general discussion of the model and the information required 
for its application is presented next. 

i. figanalMethod 

The over-all approach to the square-grid model for studying the 
effects of decontamination of fallout radiation is based on the following 
criteria: 

(a) The techniques must enable local civil defense personnel in 
the field to calculate the effects of large-scale decontamination of 
municipal areas; 

(b) The approximate techniques and procedures must not require 
computers; i.e., graphical and tabular forms should be applicable. 

As noted in Section II, the square-grid model is used to estimate 

the reduction factor D/D for the various Dianes of falliut contamination • 
o 

each at differing heights. The planes of contamination are divided into 
square-grid areas, and the detector may in principle be located at any 
arbitrary position, at height z. The size of the square-grid sections was 
chosen as 50 feet, an area which is large enough that subsequent computations 
are not unwieldy, but which is also sufficiently small to yield an acceptable 
approximation of major contaminated planes in almost any area. Consequently, 
any plane of contamination in this model is approximated by 50 foot square 
grids. An irotropic detector is assumed to be located at a height z above 
the center square grid. It is at this point detsetor that the reduction 
factor D/D q of any given contaminated plane area is computed. 

Since there will in general be several planes of contamination at 
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varying heights, contributing to the fallout gamma radiation at a given 
detector location, the analysis will functionally depend to a large extent 
on the variable z, the detector height. By varying the detector height 
above the ground plane, one can generate functions which can be used for 
the analysis of the contributions of planes of contamination at various 
heights. 

In this model, the z dependence of the detector response is treated 
separately from the other functional dependences. This is accomplished by 
separately considering the detector response from planes of contamination at 
each vertical height. The responses due to planes of contamination at all heights 
are then superimposed to obtain the total detector response or reduction factor 
D/D q at a given point. 

The model will thus work as follows: A given plane of contamination at a 
particular height is approximated by 50 foot square grid areas. The detector 
is always vertically centered on the center of the center square in the 
square-grid network. The squares are numbered as discussed in Section II. 

When the detector response due to a given square is desired, the solid angle 
subtended by the square at the detector and the effective distance from the 
square area to the detector must be computed. 

Once the solid angles are determined, the detector response due to this 
particular square may be determined by use of appropriate curves. Every square- 
grid area for the given plane of contamination may be so treated and the solid 
angles, weighting functions, and detector responses (reduction factors) of 
each of these squares may be tabulated. This procedure may be then applied to 
all other planes at another height, and so on. For the square-grid overlay, one 

then has a tabulation of the unique (in the absence of barriers) detector response 

* 

for each individual square-grid area. 

*If barriers are present, one refers to specially contructed tables and curves 
that will be discussed in paragraph 3, Section II. 
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When the square-grid overlay is appropriately scaled to detailed municipal 
maps (for example, Sanborn Maps), then the properly scaled overlays may be 
placed over the maps and the contribution to the detector due to given areas 
of contamination may be quickly assessed. This is accomplished by noting the 
square-grid numbers of the areas covering the contaminated areas of interest, 
determining the reduction factors D^/D^ for each of the individual square 
areas from the tabulated values, and summing the resultant individual 
contributions. 

D. Contents 

Section II of this appendix presents the explicit formulation of the square- 
grid model including the necessary background material and mathematical expressions 
used. Section IV illustrates the use of the model with a concrete example. An 
overlay grid is constructed and applied to a hypothetical municipal area. 
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II. EXPLICIT FOPjflJLATluN OF THE SQUARE-GRID MODEL 


A. Approach 

The notation used throughout this research memorandum is similar to that used 
in Reference C-I. The planes of contamination are divided into 5G-foot square 
sections (determined by a grid overlay) with the detector centered vertically over 
any desired grid (chosen in each instance to be the origin of a rectangular 
coordinate system). 

The detector response due to the individual square lattice grids is then 
computed using the L(X,o>) functions,* provided the effective attenuation distance 
X and the solid angle w> are known. It should be noted here, however, that the 
determination of the detector responses requires special precautions since the 
L(X,cn) curves of Spencer are all computed for a contaminated plane which is 
symmetric about the detector location. This is not the case with square-grids 
in the square lattice model. Consider the situation depicted in Figure C-l. If 
one desires to know the contribution due to the square-grid in the upper right 
hand side of the Figure, it is first necessary to compute the L(X,u>) function for 
the solid angle subtended by the entire square annular region enclosed by the 
broken lines in Figure C-l. Then an appropriate weighting must be applied to 
estimate the detector response from the individual square-grid. This weighting is 
not as simple to compute as in the case of the circular model vroare a simple propor¬ 
tion was sufficient. The weighting scheme devised fir the square-grid model it 

r it was found that theae functions could be approximated over their range of 
greatest utility by simple analytical sxpresslon* (Section Ill). Such 
expressions can be evaluated with the aid of a desk calculator and standard 
mathoaaticsl tables. 
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Figure C-l 

Concentric Squares Configurations 
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discussed in Section n-B. 


B. Formulation of the Model 

1. Assumptions 

The basic assumptions underlying the square-grid jdei are: 

(a) Each 50 foot square is contaminated uniformly; 

(b) The gamma radiation dose received by 'he detector from a shielded 
contaminated square-grid is due to both direct gaum radiation 
passing through the barrier and air-scattered gamma radiation; 

(c) The detector response from disjoint sources is equal to the sum 
of the contributions from the individual soui-es; and 

(d) It is verified that the orientation of a square ~£rid with respect 
to the detector location need not be considered. The distance 
from the center of an individual 50 foot square to the detector 
and the distance between lattice points provide the geometric 
characteristics required in the computation of detector response. 
This property is of considerable practical importance to the 
construction of simple graphs and tables for determining the 
solid angle subtended at the detector by any lattice square as a 
function of distance. 

2. geometry 

The square lattice or square-grid model is constructed in a manner 
similar to the circular model of Appendix B, The two basic differences 
between the circular model and the square-grid model lie in: 

(a) The procedures for computing the solid angles; and 

(b) The techniques for treating the position-dependent parameters 
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used to compute the individual contributions of square-grids. The 
lattice uses a unit-square annulus (Figure C-l); azimuthal sector 
contributions are used in the circular model. 

Consider the labeling scheme for the unit squares of the lattice indicated 
in Figure C-2. In this Figure, the first number of a pair (i) in a unit square 
identifies the unit-square annulus. The second number (j) specifies the 
position of the unit square in the square annulus by indexing counter clockwise. 
In addition, there is a third number associated with each unit square, h^, 
which specifies the height of a unit square above the base reference plane. 

If we let a be the width of a unit annulus, then the solid angle 
fractions subtended by the n^ square annular region is given by 

A 0 ) « t <W 2 - h n ) " T * x n-l’ y n-l’ Z "V * (C " 3) 

The solid angle fraction subtended by a unit square is given by 


“IJ * * [ V tx ij + 2 • Hi + 2 ■ 2 * V 

- q 2 y,j +f. *- » 13 > 


+ q 3 * 


’ij ' f 


, z - h ) 


- q 4 y.j -f, * - h tj ) ] 


(C-4) 


r(x,y,z - h) - £ tan 1 / -ittl—ixl 


/ MJxi_ __\ 

\<z - h) If x 2 + y 2 + (* - h ) l J 


(05) 




(C-6) 
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Figure 02 

Square-Grid Labeling Scheme 
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and the coordinates (X^, and y^) are the coordinates of the center of the 
i, j square. 

The above equations are essentially those presented in Spencer's 
Monograph 42 (Reference C-l) on pages 68 and 69. 

Differences in the L c (X,o) functions tabulated in Spencer's Monograph 

are appropriate to give the reduction factor (^“) for the square annular 

o 

regions described above. Of course, the proper solid angles must be substi¬ 
tuted in their arguments. 

In this section, barrier shielding of contributing source squares is not 
treated. This is treated separately in the next section of this report. 
Barriers require the use of the L a (X,u>) function in addition to the L c (X,oo) 
function. These, functions may be approximated by simple functions of the form 

L = Auj m exp(BX) , (07) 

where A is the intercept value and B is the slope of the curve (See 
Section III). 

In order to show specifically the scheme for computing the detector 
response, Figure 03 is used. 

In terms of this notation and the discussion of this section, it follows 
D 

that the ^ of this typical square uniform annulus region is given by 

£ * [MW - Wi’Vi'] L(X) • < c - 8 > 

Since there are 8n unit squares in the i> C ^ square annulus (n £ 1), the 
contribution of a typical unit square of the square annular region may be 


012 
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written in the form 



JLl 

D 


i% [ L c (X n,j,V- L c (, .-l,J,Vl>] 


L(X) 


(C-9) 


8n 8n 


(C-10) 


The factor o> n j/^ and u> n * £ co n ” o> n _^ is the weighting which 

corrects for the effective position of the different unit squares of a given 
square annular region. 


It is easily seen that the contribution from each square of a given 
square annulus n is weighted according to the solid angle fraction it 
subtends. The total geometry reduction factor for the whole plane is thus 
given by 


$- ■ ? ? 


? L (X . a) ) - L (X . . a) ,) 

j 8n« n c' n,j, n 7 c v n-l,j, n-1' 


L(X) (C-ll) 


D - ' I ? Z 

o n j 


£ [ l c<Vj,V - VVi.j.vi*] 


L(X) . (C-12) 


In order to illustrate the manner in which the square-grid model is 
applied, an overlay was constructed of 50 foot square-grid areas (Joe Figure 
C-2). In order to us the overlay for decontamination analysis, the reduction 
factor, D/D q , for each of the 50 foot squares is required. The overlay scaled 
to appropriate local maps could then be quickly and effectively utilized in 
accessing the effects of large scale decontamination in a particular area of 
interest. In order to obtain the D/D q values, the solid angle of each 
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Individual 50 foot square must first be determined. An appropriate weight is 
then obtained by computing the total solid angle of the annular region in 
which a particular 50 foot square is located; the solid angle fraction of the 
individual squares may then be determined. 

The individual calculations of the solid angles are quite, lengthy. In 
keeping with the spirit of the model, an approximate but highly accurate 
procedure for computing the solid angles as a function of horizontal distance 
is derived for various detector heights. Computations of the solid angles 
of the 50 foot square-grids were made for appropriately spaced on-axis center 
distances. (The on-axis case has desirable symmetry properties which simplify 
the calculations.) To illustrate this construction, » defector height of 
100 feet was used (z = 100 feet). Calculations of the solid angles subtended 
at the detector by the 50 foot squares were made for each of the successive 
on-axis square-grids. Calculations for several off-axis lattice squares were 
performed (e.g., x a y = 150 feet). It was found that the solid angles of the 
50 foot squares depended (to a high degree of accuracy) only on the absolute 
distance of the square lattice from the detector. With z given, the 
dependence is a function of horizontal distance. Consequently, the curve of 
Figure C-4 was drawn showing the approximate solid angle of any square lattice 
for three detector heights. The solid angles read from Figure C-4 are given 
in Table C-I for z ■ 100 feet. The solid angle fractions used lor 

weighting tne L(X,u) functions are also given. These solid angle fractions 
are computed by use of other calculations summarized in Table C-Ii. "he 
reduction factors D/D q (no barriers) for a particular square annular region, 
n, at a horizontal distance (X) from the plane of contamination is then given 
by 

I? ' «» [ l c«-n> ' L c (X ’ a-l>] 
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Solid Angle Fraction 


3 



(feet) 


(Horizontal Distance from Detector) 



Figure C-4 

Solid Angie Fraction of SO Foot Square* aa a Function of 
Horizontal Distance froe Detector 
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TABLE C-I 


Solid Angles and Solid Angle Fractions for 
the 50 Foot Square.-., (z a 100 Feet) 


HIBfPSSH 

MM| 


mi 

HP! 

^L(7.6) a. c (7.6, % ) 

0.0 

0 

_ 0 

mu 

1.00 

.00632 

i,i 

50 

0 

.02775 

,1397 

.00620 

1,2 

50 

50 

.0212 

.1077 

.00478 

2,1 

100 

mm 

.01416 

.0774 

.00399 

2,2 

100 

MBm 

.0122 

.0667 

.00346 

2,3 

100 


.0079 

.0432 

.00223 


150 

mm 

.00689 

.0547 

.00224 


150 

||BjH 

.00620 

.0493 

.00202 


150 


.0046 

.0366 

.00150 

WmZM 

150 


.0032 

.0254 

.00104 


200 

0 

.0036 

.0421 

.00148 


200 

50 

.0034 

.0397 

.00139 


200 

100 

.0029 

.0339 

.00118 


200 

150 

.0022 

.0257 

.00090 

mam 

200 

200 

.0016 

.0187 

.00066 

77^ 

250 

0 

.0022 

.0366 

.00107 


250 

50 

.0021 

.0350 

.00103 

wr;- « 

250 

100 

.0018 

.030 

.00087 


250 

150 

. 0014 

.0233 

.00068 


250 

200 

.0010 

.0167 

.00049 

■ 


250 

.00079 

.010 

.00029 

6,1 

300 

0 

.00127 

.0288 

.00051 

6,2 

300 

50 

.00120 

.0272 

.00048 

6,3 

300 

100 

.00105 

.0238 

.00042 

6,4 

300 

150 

.0GQ8 

.0182 

.000338 

6,5 

300 

200 

.00063 

.0143 

.000255 

6,6 

300 

2 50 

.00045 

.0102 

.000180 

6,7 

300 

300 

.00040 

.0091 

.000157 
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TABLE C-II 

Solid Angles Subtended at the Detector (0.0.100 1 3 
bv the Square Annular Regions Centered at the Origin 



w n 

Atat=ta> -tat , 

n n-1 

n “ 0 

0.0375 

0.0375 

l 

0.2344 

0.1970 

2 

0.4i;5 

0.1830 

3 

0.5436 

0.1261 

4 

0.6292 

0.0856 

5 

0.6892 

0.06C1 

6 

0.7333 

0.0441 


These reduction factors for various annular regious are given for z * 100' in 
Table C-I1I. With this information, the contribut:(and hence the effect of 
decontamination) of each of the square-grid areas shown in the overlay diagram 
(Figure C-2) can be tabulated. Consequently, the overlay can then be placed on 
a scaled city map and the direct contribution of any area can be <-?siiy 
determined from the tabulated contributions. A table of contributions can be 
prepared for any desired detector height, z. For large-scale considerations, 
the 50 foot square areas are small enough to be linearly divided for partial 
contributions with little error. 


TABLE C-II 


Sam, .imlsjL3ub&ssi$s£. 


fie j*te ctoy , ( 0J) JLgo.', 1 


hauare Annular Region 

D /D - L(7.6)X. (7-6,.v ) 

a. Q . ... . v.ft . 

n * o 

0.00b 3 

l 

0.0**44 

2 

O.C315 

3 

O 0409 

4 

0.0351 

5 

0. 

b 

0.v)17b 
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3. Barriers 


An approximate treatment of barrier attenuation and scattering is mm 
considered. Since the model is constrained to stay within the framework of 
the results of Spencer (Reference C-l), a formulation is developed in terms of 
Spencer's L c (X,w) and L a (X,w) functions (Reference C-l). Here, the detector 
is assumed to be located in an open area. Howevet, this assumption ^ocs not 
restrict the application of the model lo regions in which the detector is in 
an open area. If the detector is in a structure, it will be necessary to 
analyze the shielding by that structure. Procedures both approximate and 
more detailed are available for this analysis (Reference C-2). 

The purpose of this section will not be to enter into the treatment of 
barrier problems in general. Rather, the purpose, is to characterize the 
treatment of barriers in a manner that is consistent with the approximate 
character of the over-all model and which is simple enough to be practical 
for field applications. The formal treatment of barriers in the spirit of 
the above discussion is first considered. 

Consider the contribution D n j/D Q of the unit square of the n^ square 
annular region to the detector located at the point (x=0, y=0, z) in the 
absence of any barriers: 

0 f L (X . ) - l (X , a . J I L(X) (C-13) 

nj o 4i'> n I c' air, n c' air n-1 j 

vher; X . is the effective mass thickness of air between the detector and 
air 

the source plane. If barrier* arc present which intersect the solid angle 
cone subtended at the detector by the tyj** 1 square, the direct and the 
scattered components of the radiation arc affected. With intervening 
barriers present in the n, n-1, .... square annular regions, we have, 

in principle, a total s«ss thickness f X^ ^ defined vim the respect to the 
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n,g square and the given detector location, i.e., 

t X n,g ~ a n,n,i X n,i + ^ a n,n-l,i X n-l,i + ... + ^ a n,n-m,i X n-m,i + 

1 f '’-‘n.J.i) X .ir 



fraction or the solid angle (defined with respect to intersected by the 

barrier at the n,j t ^ square. The last term of the abpvd^'expression treats 
approximately the unobstructed fraction of the solid angle which is filled with 
air It should be noted that in practice there would probably never be more 
tha.'. three terms, and, in general, fewer, in these sums over the mass thickness 
parameters cf the barriers located in any square annular region. 


derived. With u barrier of effective mass .hickness, ^ ^ symmetrically 
sunounding the detector between the detector and the n C ^ annular region. 

V D o- I V„,j> [VA.W • VA, j-vi>] ■ < c - l4 > 


The fractional part (D^/D^) is not the desired contri¬ 

bution because this fraction: 


(b) does not include the scattering arou. the barriers of radiation 


originating in 


n » J 


020 




usual form to compute this contribution because they include a large component 


of scattering into u> n ^ from the infinite plane beyond u^ j and no contribution 
from the shielded area which scatters and arrives at the detector outside of 


a n,j‘ 

However, by defining a function in terms of the L a (X,co) in a manner 
analogous to that in which L c (X,u>) is lefined in terms of L(X), we obtain a 
function approximating the contribution from a given area. 

L(—&-) 

fit L »c. (X ’ m> ■ L a <X -> ' • 


This function eliminates the contribution of scattering into « n ^ fron- 

the infinite portion of the plane not included in u> n ^. Differences of the 

L fic , (i.e., AL flc ) include a scattering contribution into the solid angle Au.^ 

from that portion of the source plane corresponding to the area defined by 

the solid angle a ,. The AL functions also include the direct contribution 
n-l ac 

from the plane into A to . However, the AL does not include a contribution 
v n ac 

to the detector by radiation originating in “n-i which is scattered out 
v. u> n * co n _^ and subsequently scattered back intc the detector. 

These are approximately compensating contributions and we will use this 
function to compute the approximate contribution from barrier obstructed 
source squares. x Clearly, for situations where most of the radiation from 
any given source square is unobstructed by barriers, this approximation will 
be adequate. When most of the source planes are obstructed by very large 
barriers, then (3.2) would be a better approximation. Thus, we cake for the 
approximate contribution of the n,j th square source area with effective mass 


’fore exact sxpreasiona could be developed for this contribution; however, it 
waa felt -hat Chla effort would not be Justified until after the potentiali¬ 
ties of a square-grid point aourc* model sd been explored. 
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thickness, t X R y the expression 

* Sf L <t X „,J> ^.cVn,. “n> 

' ZSf L <A,J> [ L .c<A,J-“n> ‘ L a, ( t X n,j'“ i n-l ) ] • 


(C-15) 


We will show below that the values of (C-15) may be presented 
graphically for a large range of X. This represents a conservative approxi¬ 
mation when most of the source squares are not obstructed by barriers. 

The detector response due to a given contaminated 50 foot square area is 
given by the function 

f(Q)L(X)AL ac 

where AL is defined by the relation 

' L .c (1 W - 


with 

L .c * Trn [w>V x i’“) - ■ 

The quantity f(u)) is the Appropriate weighting for a given square grid area 
and ia given by a/ds. The function L(X)AL ac is highly dependent on the 
barrier mass thickness X. However, for a given annular square, (i,j), 
ft-)L(X)Ai. #c is dependent on the solid angle of the Individual square only 
through the quantity f(x). The quantity f(u) varies with the position cf the 
50 foot lattice in a given annular square i. Consequently, the quantity 
L(X)AL at may be computed as a function of X for any given square annulus. 
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Then, through the f(co) function, the reduction factor of each of the individual 
50 foot squares penetrating a barrier of any thickness S may be easily 
obtained. The function L(X)AL ac as a function of X was computed for rings 
0 through 6, and the results tabulated in Table C-IV. Figure 05 also 
presents these data graphically. It is noted that the variation of 
L(X)AL ac as a function of X is for all practical purposes a single, 
exponential. (The variations from the straight line on the semi-logarithmic 
plot are well within the accuracy witn which Spencer's curves L(X) and L a (X,to) 
can be read. 

Table OIV gives the values of f(ce) - co/Acjo for each of the 50 foot square- 
grids through ring six. Although these weighting functions for any given 
annular sector, i, vary from one square-grid to another, an average value 
of f(co) for each ring can be obtained chrough the function Aco8i, where 8i is 
the number of 50 foot square-grids in ring number i. The accuracy of using 
the average solid angle for every square-grid in an annular ring is certainly 
well within the general accurcy obtainable for the reduction factor in any 
problem in which complex barriers are present. Thus, we assign this average 
value to each 50 foot square-grid in a particular ring. With this assignment, 
the reduction factor of any square-grid through a barrier of arbitrary thickness 
may be computed and presented ir. graphical form. Table C-V presents the 
necessary data from which Figure 06 is drawn. Figure 06 thus represents 
the reduction factor D/D q for the individual square-grid areas as a function 
of the barrier mass thickness X. 

Thus, the ease of readily superimposing the- values of D/D^ for the 
various square-grid areas for the case of air with no barrier is now carried 
completely over to the case with barrier. One needs three Ingredients: the 
barrier mass thickness X; the square-grid overlay; and the graphical repre¬ 
sentation of the function* ffor the various rings, for example, 
Figures 05 and C-b. 
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TABLE C-IV 
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Figure 0 5 

Variation or L(X)AL^ c vs. X 
For Ring i - J 
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Figure 06 

Reduction factor D/D 0 for the 
Individual Square Grid Area* 
as a Function of the Barrier 
Mass Thickness X for Ring 
1-3. 
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The expression (C-14) which gives the total effective mass thickness from 
the detector through the intervening barriers to the n.j^ square-grid area 
is in -,c?neral quite difficult to evaluate. However, one can argue on physical 
grounds, that, for the regions of interest In the particular applications of 
this model, the computation of the effective mass thickness can be greatly 
simplified. Indeed it is possible to show theoretically that the expression 

(C-14) is approximately related in a rather simple manner to the effective 

th 

line distance from the detector to the n,j square area for large distances 
and barrier detail in keeping with the approximate character of this model. 
Then, the total effective mass thickness £ X n ^ can be practically estimated 
as follows: 

(a) Compute the line distance from the detector to the 
center of the n,j C ^ square radiation source area; 

(b) From considerations of the barriers (i.e., wall, floor, 
roof, etc., data) collected from earlier building 
surveys, compute the effective distance through which 
the radiation travels in the barriers; i.e., simply 
the thickness of all the walls through which the 
radiation travels (note: the obliqueness cancels 

when the area is projected to the vertical). 

(c) Add the effective thickness for each type of barrier 
material separately. Then use the formula for X 

in conjunction with the material's mass thickness 
table (see Reference C-l, p. 15) for each material, 
then sum X's for all the materials to obtain 
*.he total barrier mass thickness. 
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III. APPROXIMATE ANALYTICAL FIT FOR L^X.oq) AND L c (X f tu) FUNCTIONS 


It was found during the course of developing the circular model (Appendix B) 
that an analytical fit over a rather extensive range of values for both the 
L fl (X, oj) and L c (X,oo) curves found on pages 41 and 42, respectively, of Reference 
C-l is possible by employing an expression of the general form 

D/D q * L(X,a)) = a 8 e** , (C-15) 

where s is the slope and A the Intercept for selected values of 
X and to. The slope in turn is obtained by 

s = In L(X,cd 2 ) - In Ltf.u.j.) , (016) 

while the intercept is calculated by solving Equation C-15 for A as 
follows: 

A = la itULml - 5 la-m 

x 

The best analytical fit for the L c (X,il) curves was achieved by obtaining 
the L c (X, ) values given on page 41 of Reference C-2 corresponding to ^ and 
u '2 values of 0.05 and 0.2 respectively and solving for the intercept at , • 0.2 
for each X value selected. The same procedure was followed for the L # (X,..) 
curves depicted on page 42 of Reference 07. The slope in this case was 
calculated for ■ 0.01 and - 0.1 while the Intercept was obtained at 

• 0.1. The numerical results of these calculations are presented in Table 
C-VI, -nd the resultant curves are shown in Figures 07 and 08. 


It is apparent when comparing th*. L(X, . ) curve* developed here with those 











































































appearing in Reference C-2 that the accuracy of the L(X,od) values Deteriorate 
rapidly due to the lack of convergence as to + 0. This is particularly true 
for those values of us > 0.3 and cd >0.5, associated with the L fl (X,a)) and 
L c (X,aO functions respectively. Despite the accuracy limitations inherent 
with this analytical approximation of L(X,cu) functions, there is no evidence 
that the approach would not produce the desired results when appropriately 
used in the square-grid model. It should be noted that there are no 
restrictions associated with this model that would prohibit the exclusive use 
of the exact values for L a (X,to) and L c (X,o>) presented in Reference C-l. 


TABLE VI 

galcuiateQalues_cf_s_and ._A „ ,. for 

Selected Effective Mass Thickness Values 


L a (X,a>) 

L c (X,u>) | 

X 

s 

A 

x 

5 

A 

1 

1.125 

-0.725 

2.26 

1.115 

-0.396 

5 

1.055 

-0.0950 

7.53 

1.137 

-0.0531 

24 

1.055 

-0.0952 

22.6 

1.042 

-0.002,55 

72 

0.975 

+0.00445 

75.3 

0.916 

40.00726 

144 

0.957 

+0.004355 




216 

0.954 

+0.00361 




288 

0.930 

+0.00302 
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IV. AN EXAMPLE 


Consider the situation indicated in Figure 09 below. We wish to explicitly 

demonstra i the method of calculation of the D ./D for several of the contri- 

n, j o 

buting squares of the grid. Starting with the application of the overlay, 
discussed in Section II-B (see Figure C-3), to the quadrant of the idealized 
scaled map shown in Figure OS for all of the contributions D r j/D q of the squares 
can be obtained directly from the appropriate overlay as discussed earlier. Of 
course, the values will only be correct when the radiating source square contri¬ 
buting to the detector is not appreciably obstructed by barriers. In particular, 
for the contribution of the 2,3 square (for a detector located 100' above the 


L (7.6,0.234)] = 0.00218 is immediately read off of the overlay. 


We now assume that an equivalent thickness of 2' of concrete for the building 
blocking the contribution of square (3,7) has been determined. Also, we assume 
that an equivalent thickness of 1' of concrete for the same structure blocking 
the contribution of its roof, square 2,5 has been determined. Then the effective 
mass thicknesses of these barriers are given by (see Reference C-l, p. 15) 


X 3 ? = 2(z/A)pA ■ 2(0.5)(144)(2') * 288, 

X2 5 • 2(0.5)(144)(1') - 144 

The approximate D n j(X)/D q for the source squares obstructed by barriers cannot 
be read from the overlay because the L tc (X,uj) function. 1 leading to this contri¬ 
bution are functions of the various effective barrier thicknesses X present. 

An important but easy-to-apply approach has been discussed for obtaining the 

terms D ,(X)/D in Section II, which employs a table and/or graph giving the 
n, j o 
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Figure 09 

Idealized Scaled Hap (Co be uzed in conjunction 
with Overlay - tee Figure 03 end Table OI) 
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values of the obstructed contribution of each source square D /D as a function 

n,j 0 

of the effective mass thickness X (see Table C-V and Figure C-5). Using the 
graph or Table C-V to evaluate D. „(X = 288)/D for the example of this section, 


formal expressions representing the approximate barrier obstructed contributions 
of the source squares 3,7 and 2,3 (Detector location (0,0,100')) are given by 


3,7 / o 


L(288) L ac (288,o) 3 ) - L &c U88,c 


D 2,5 /D o- ^ 


L(144) L (144,^ - L (144,0),] 


Using the information on the overlay relating to these solid angles we can 


D 2 5 /D q » (0.077) L(144) L #c (144,0.418) - 1^(144 


,0.234) 


Then using Spencer's tables (Reference C-l, pp. 37,42) we find: 


- (0.055) | (0.00024)(0.067)} 


jjk* - (0.077) { (0.0063)(0.349)| 


C-35 


and 

D 2 5 /D o “ 1,54 x 10 " 4 

Of course, these results could also have been obtained from a table of the 
function L ac as a function of to with X as a parameter. 
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V. SUMMARY AND CONCLUSIONS 


It has been shown that a square-grid overlay technique may be used to analyze 
the radiation dose contribution to a detector from individual planes of 
contamination. Once these contributions are known, the effectiveness of 
decontaminating any plane may easily be assessed. Data required for appli¬ 
cation of tne square-grid model are; 

(a) A grid overlay with each square labeled in a coordinate 
scheme (a 50 foot grid size is proposed with each square 
labeled (i,j) where i is the square aiomlus containing 
the square (i,j) and j is the position of the square 
(i,j) in the annulus. The detector is vertically above 
the center of square (0,0)). 

(b) A map of the area under analysis which is of the same 
scale as the grid overlay. 

(c) Data on the effective mass thicknesses of structures in 
the area. 

(d) A set of tables giving the contribution from each square 
in the grid to the detector located a given distance 
above the plane of the grid. The tables are constructed 
as a function of the detector height above the plane. 

(e) A set of curves giving the effects of barrier shielding (as 
a function of X) on the detector response for each square 
in the grid. 

The required tables and curves are derived for a detector height of 100 feet 
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•r 


and a square-grid length of, 5C feet. It is shown for this case that the grid may be 
approximated by a function of the distuice of the center of the square from the 
center of the (0,0) square. Changing the orientation of the square with respect 
to the detector produces only a very small perturbation in the value of the solid 
angle. 

It is concluded that when detailed maps of a standard scale are available 
for an area (such as Sanborn maps) a square-grid overlay technique can offer a 
practical tool for analysis of decontamination effectiveness. It is recommended 
that data be derived for a useful range of detector heights (say 3', 10', 30', as 
well as the present 100') and that the method be applied to some real situations. 

Further work in barrier effects is required. A more refined approach to the 
analysis of barrier effects should be .sought. 
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GLOSSARY OF TERMS 


Ap 

%Ap 

Centered area 


»„<V |h t -z|, 
w !h i- 2! - r) 
VV 

B v <V 

VlV'l’ r t> 
c d < l h i" , l' V 

e 

h i 

L(X) 

r i 


“ area in square feet of the l C ^ contaminated area 
- fraction of the vertical barrier voided by apertures 

■ percent of the vertical barrier voided by apertures 

■ an area of contamination such that the line from the 
detector to the center point of the area is perpendicu¬ 
lar to the plane of the area 

■ horizontal barrier shielding factor for centered areas 

■ horizontal barrier shielding factor for off-center area 

* vertical barrier shielding factor (no horizontal 
barrier between source and vertical barrier) 

■ vertical barrier shielding factor with horizontal 
barr -s between the source and vertical barrier 

= doee cor ribution per unit area 
« doee contribution from centered area 

* eccentricity ratio, length/width 

*= height of the i C ^ area above reference plane 
«* total dose from an infinite plane source of fallout 

■ horizontal distance from the detector to the center of 
the i th area 

■ solid angle fraction 

* horizontal-shield equivalent mass thickness in psf 

■ vertical-shield equivalent ma^s thickness in pcf 

■ height of detector above reference plane 
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Appendix D 


A Point-Source Model for Approximating 
Gamma Ray Intensity at a Single Detector Location 

I. INTRODUCTION 


A. Purpose 

This appendix presents a procedure for estimating the relative radiation dose 
contributions at a detector from multiple planes of contamination. The two 
objectives in developing this model are: 

1. to provide accurate approximations to the relative dose contributions 
by many different planes of fallout contamination, and 

2. to reduce to a minimum the charts and tables, maps, equipment, and time 
required to calculate the contributions to a given detector location. 

B. Background 

Two other models with similar objectives have been developed under this 
contract. Conceptually, the "circular model" (Appendix B) is simple, and it 
is easy to use, but accurate approximations are difficult to make when the con¬ 
taminated planes under analysts are not segments of circular annuli. The "square 
grid model" (Appendix C) can be used to approximate accurately the more usual 
geometric arrangement of rectangular planes, but application requires a suitable 
scale map for use with overlays of che grid. This point-source model presents 
a scheme for approximating contaminated planes by the appropriately located 
point-sources of radiation. The point-source model is formulated such that 
accuracy limits are known. Scale maps are not required, but only certain dimen¬ 
sions; only a small number of concise charts are required for its application. 

C. The Model 

In this model the detector and source point geometry are presented in cylin- 






drical geometry. Since the detector is isotropic, the geometry reduces to two 
distances. These are |h-zj, the vertical separation of the detector from the 
source plane, and r , the horizontal separation of the detector axis from the 
point representing a given contaminated area. Here h is the height of the source 
plane, and z is the height of the detector, both above a reference plane. Figure D-l 
illustrates this geometry. 



Figure D-l 
Geometry of Model 

The shape of the contaminated areas to be represented by a point is of con¬ 
cern. It is shown that for the area centered vertically below (or above) the 
detector, a circle or square of a given area subtends essentially the tame solid 
angle and thus yields the same relative dose. Off-center areas are all assumed to 
be squares. Off-center areas which are not square should be approximated as closely 
as possible with squares, unless the maximum dimension meets the criterion of 
Figure D-ll. 
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It has been shown (Appendix C) that the orientation of a square area within 
a source plane does not appreciably affect the solid angle subtended at a detector. 
Thus, to a good approximation the relative contribution by a contaminatec area to 
a detector is a function only of the vertical height |h-z| separating source and 
detector, the horizontal distance r from the center of the area to the vertical 
detector axis, the size of the area, and the intervening shielding. 

In this paper the relative dose rate, or reduction factor delivered to the 
detector from a given area is developed. This relative dose rate for area i is 
given by the ratio Dj/D q where is the dose rate contributed to the detector by 
plane i and D q is the total relative dose rate three feet above an infinite, 
uniform plane of contamination. The reduction factor RF is then given by 

RF = 2 D /D 

all planes 1 ° 

The relative contribution from all areas beyond those analyzed individually may 
be obtained by subtracting the "centered area" relative contribution, Figure D-3, 
from the infinite plane contribution L(|h-z|) given in the right margin of 
Figure D-3. 

A major consideration in the development of this model was the establishment 
of the maximum area of a square which could be accurately represented by a point- 
source at the center of the square. A criterion was devised and maximum areas 
were calculated as a function of horizontal separation for several vertical 
separation distances (|h-z| » 3', 10', 20', 50', and 100'). The criterion 
establishes that for a square area less than or equal to the maximum area, the 
relative contribution calculated on an area basis from Spencer's curves (Reference 
D-l) are within 107. of the contribution calculated from the equation 

C * A i’1 h i"*l ,r i * C 1 ^ , i’*l ,r i* A i (D-l) 
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where (|h^-z|,r^) is the contribution per unit area as a function of |h^-z| and 
r^ and where is the contaminated area in square feet. 

The factors for the shielding effects of barriers between the source and the 
detector are derived from Spencer's curves and other established sources (Ref¬ 
erences D-3, D-6, b. D-5). They are presented in the standard form of multiplicitive 
adjustments of the unshielded dose from a source. 



II. DEVELOPMENT OF THE POINT-SOURCE MODEL 


The analysis of the effects of decontamination of multiple contaminated 
planes surrounding a detector usually presents unique requirements. This model 
seeks to provide a concise, rapid method of calculating accurately the fraction of 
the total fallout dose rate at a detector contributed by each contaminated area. 

In decontamination analyses, the detector positions of interest are usually unpro¬ 
tected or are in low protection factor structures. This is because decontamination 
is most effective in hastening the resumption of operations in poorly protected 
areas (Reference D-2). Shielding is treated in an approximate manner consistent 
with its relative importance in the model and with the requirement of conciseness. 
Techniques of the "Equivalent Building Method" (Reference D-5) and the "Engineering 
Manual" (Reference D-3) are combined. 

A. Areas Centered Opposite Detector 

The area immediately below or above a detector is always of interest in 
determining radiation doses. This area may be a major contributor to the dose in 
the absence of significant shielding. In any case, the difference between a 
centered area dose and the infinite plane dose (assuming uniform, unshielded 
contributions) is required as the contribution from beyond the areas which are 
analyzed specifically. 

1. Centered'Circles and Squares 

In this model we wish to approximate the dose from centered areas by a 
function of only area and detector height. If we are tc eliminate shape as 
a parameter, the most obvious concern is whether or not centered circles and 
squares of equal areas contribute equally to the detector. Assuming that 
the contribution from a .entered circle or square is a function only of 
letector height and o f the solid angle subtended by the area, the desired 
approximation of an equivalence criteria between squares and circles is that 




equal solid angles are subtended by equal areas at equal relative detector 

heights. Comparisons of solid angle fractions subtended are given in Figure 

D-2, It is seen from this figure that for the range of areas and detector 

heights used in this model, tho solid angle fraction subtended by centered 

circles and squares of equal area are equal within the limits of accuracy which 

can be used in reading dose rate contributions from graphs. Thus, the dose 

rate contribution from centered circular areas is taken as being equal to the 

dose rate contribution from centered squares of equal areas, as is the 

standard practice (References D-l and D-3). 

In Figure D-3 the relative dose rate contributions from a contaminated 

circular or square area centered above or below the detector is presented as 

a function of area for several detector heights |h^-z|. These values were 

derived in exactly the same manner as were Spencer's L (X,u) functions where 

c 

X is the equivalent mass thickness between the source and detector and u 


is the solid angle fraction subtended by the source (Reference D-l). That 


Ih^-zj feet above an infinite plane of fallout contamination of such intensity 


that L(3') * 1, then 


C d(l h t - z I,A t ) - L(|h. 


i-*l> ' L (fV 


^|h i -z | 2 + r i 2 « | h^ -z | /1-u^ m 


Thus, c d (|h^-e(, A^) corresponds to Spencer's L(X)L c (X,u) function. The 
L( |h^-*|) values given in the margin of Figure D-3 are the relative dose ri 
from an infinite area of fallout contamination, A » », at a detector at 
height |h t -*j . 
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Figure D-2 

Comparison of Solid Angles for Centered Areas 














Centered Under Detector 


















2. Centered Rectangles 


Areas encountered in practical situations are roost often rectangular. 

We shall thus investigate the maximum eccentricity of a centered rectangular area 
such that the dose rate contribution is within a giver margin of the dose 
rate from a centered square (or circle) of the same arec. 

A centered rectangular area (i,e., L/w + 1) will always contribute less 
dose to the detector than a centered square of the same area. Thus by approx¬ 
imating a centered rectangle by a centered square of the same area, one will 
overestimate the dose contribution from the area and consequently will over¬ 
estimate the effectiveness of decontaminating the area. This means, of course, 
that the dose rate overestimates should be kept to a small percent of the doss 
rate contribution from the area. This overestimate of the relative dose rate 
of decontaminating the centered area will mean that thr reduction of decon¬ 
taminating other areas will be underestimated. We shall establish bounds on 
the error introduced sufficiently stringent to insure acceptable accuracy 
in the application of the modal. 

Considering the precision of other factors in the total decontamination 
model, an error limit of ten percent seemed acceptable. The criterion estab¬ 
lished below assures that the error in unshielded dose contribution introduced 
by this approximation is less than ten percent. 

Consider the two geometries of Figure D-4. 


This 

Configuration 

Yields 

Cy,i A~ Ai ,) 




Figure r>-4 

Approximately Equivalent Areas - Circular Segments 


D-9 



We see that 



AA - IM 
A? - M/T 

-r -V(s + "OH) 2 

where 

A • area of square "center" of rectangle 
A + AA ■ urea of rectangle 
T * width of rectangle 
T + A0 * length of rectangle 
Thus, the eccentricity e of the rectangle is 

e ‘ ’ l+S A • (B-3) 

We have shown (Figure P-2) that the dose from centered squares and circles of 
equal areas are identical within the accuracy of the graphical data. There¬ 
fore, let denote the area A in a centered circular or square configu¬ 
ration, and let A^ denote the area A in a centered rectangular configuration 
Thus, the dose rate from the area on the left in figure D-4 is C d ( [A + M] R ) 
2 

Letting A£ ■ n r' , the dose rate from the area on the right is 

c d ( [* + * W * W + the (W - W’ • (n - 4> 

c 

Since the areas are equal, we make the approximation 

c d <( A + ^]r> i <f A + M) cs ) • (P-5) 

Now consider the way in which the off-center M is approximated in Equation 
9-5. Figure D-5 depicts this approximation, and it is clear that the circular 
segment approximation of Equation D-5 uses areas farther from the detector and 
thus gives a low approximation to the dose rate from (A + M) R . 
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Figure D-5 

Approximately Equivalent Areas - Circular Segment 


Figure D-6 depicts the way in which areas are equated in the approximation 
C d ([A + C d ([A + M] c ). (D-6) 



Figure D-6 

Approximately Equivalent Areas - Centered Square 
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In this case it is clear that the centered square approximation, Equation P-6, 
uses areas closer to the detector and thus gives a high approximation to the 
dose rate from (A + M) R • Therefore, 

C d ([A + M] c ) > C d ([A +M] r ) > C d ([A +M] ct; ) 
and 

C d ([A + AA] c ) - C d ([A + M] r ) < C ([A + M] c ) - C d ([A + M] c8 ) • 

C d ([A + AA] c ) Cj ([A + M] c ) 

The dose from the rectangular area C d ([A +M] r ) i3 then within ten percent 

of the estimate C, ([A + M] ) if 
d c 

C d ( I A - C J » A *^ CS > < 0.1 . (D~?) 

c d (U + M] C ) 

Satisfaction of Equation D-7 is thus taken as the criterion for the maximum 
eccentricity of a centered rectangle of area A + M to be approximated by 
a centered square of area A +M by using the dose curves in Figure D-3. 

The locus of maximum allowable eccentricity e from Equation D-3, as 
a function of total area, A + £ji, is derived as described belcw for the 
detector heights of 3, 10, 20, 50, and 100 feet. To obtain a point on the 
locus for a given height, the following steps are taken: 

(a) an area A is selected for the center square 

(b) incremental. AA's are added to A and the ratio of Equation D-7 
is computed for each M 

(c) the maximum M satisfying. Equation D-7 is found graphically 

(d) the maximum allowable eccentricity • ^ is calculated by 

Equation 0-3 for the area A + M. 

The values thus obtained are shown as functions of area A + , A A 

in 1'igurc D-7. The form of these curve# Is anticipated by noting that 

Um | C d ((A + ,A) C ) - C d <[A + - 0 4 

A-* « 
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Area (fc ) 


To Keep Doee Rate Error < 10% 



Figure D-? 

Maximum Eccentricity Ratioe For Centered Rectangular Areae 
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since as the center square area increases, the dose rate from A approaches 
the dose rate from the entire plane. And 


Lim | C d ([A + M] c ) - C d <[/ tM 1 g ) | - 0 

[a + m]-»o 


since a zero area contributes zero >se te regardless of the "eccentricity*'. 


B. Off-center Areas 

As implied by the title "Point-Sou Me :1", we wish to establish a scheme 
to represent areas of contamination by i int : ources. There was no requirement 
to introduce point sources for centered areas in the previous section. For off- 
center areas we shall establish criteria by which off-center areas may be represented 
by an equivalent point source at the center of tie area. Actually we will not use 
point source data. To approximate the dose contribution from the area, we will 
multiply the area by the contribution per unit area at the center of the area. 

This approach requires two sets of data for its application. First, the dose rate 
contribution per unit area as a function of horizontal distance from the detector 
is required for each detector height. Next, a set of criteria for the maximum 
area which can be treated in this manner must be established. 

1. Dose Rate Contribution Per Unit Area 

The dose rate contribution per unit area as a function of horizontal 
distance from the detector was calculated for each detector height using 
the equation 


Cj (|h-z|,r) 


C d (|h-z! ,r + ^L) - C d (|h-z|,r-~) 

—— 


(D-8) 


where 

jh-z| 

r 

C { (jh-z j,r) 


vertical separation of detector and source area 
horizontal separation of detector and center of source area 
contribution per unit area at center of source 
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(3j) a uo?a*a»<l»s i»3uoajaoH 
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Do** Rat* Contribution Per Unit Area 





















(|h-z|,r) - contribution from a centered area of radius r as 

defined in Equation D-2, (the change of argument from 
A to r is only for clarity in Equation D-8) 

2irrAr = area between circles of radius r-“ and r-t~ . 

The function C^(|h-z|,r) is presented in Figure D-8 for the selected |h-z| 
values of 3, 10, 20, 50, and 100 feet. 

The differencing in Equation D-8 to obtain these values must be accomplished 
before the maximum allowable Ar is established. To assure the validity of 
these values, the increments are chosen well below values which previous work 
(Appendix C) and judgement indicated would be the maximum values of Ar. 

These judgements are confirmed by the maximum values of Ar established below. 

2. Maximum Area to be Represented by a Point 

Since there was obviously some limit to the validity of the point source 
approximation, the following limiting arguments were devised. Consider the 
convex function shown in Figure D-9. 



Figure D-9 

Approximations on a Convex Function 
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which means that the error introduced in approximating the integral by f(r)Ar 
is less than or equal to the difference between f(r) and the average of 
f(r-|£) and f(rt^) . 

The dose per unit area functions, (|h-z|,r) are convex over most of 

their range. The log-log scale of Figure D-8 obscures this fact, but the 
general form of each of these curves is like that in Figure D-10 when 
shown on linear scales. The value r Q , such that r > r Q , implies f(r) is 
convex, is less than ten feet for all of the values of jh-z| < 100 feet. The 
inequality (D-9) thus is used in establishing the criterion for maximum area 
to be represented by a point. Only the convex portions of the (|h-*l,r) 
functions are used. 
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Horizontal Separation, r 
Figure D-10 

Form of Cj (|h-z|,r) Functions 


As stated earlier, in accordance with our objective of making a slightly 
conservative estimate of decontamination effectiveness, we establish the 
criterion 


C; (|h-z|,r-~) + (^ (|h-z | ,r-*~) 


It has been shown that the dose contribution from an off-center square 
area is insensitive to the orientation. Thus we set the maximum area to 
be represented by a point located at its center to be a square of area 
2 

A (Ar ) , where Ar is the maximum Ar satisfying the condition of 
max max max 

Equation D-10, and where r is the horizontal distance to the center of the 

square area. For areas which are not square, the criterion should bo that 

the maximum dimension is less than or equal to A* The locus of maximum 

max 
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square area to be represented by a point at the center is given in Figure D-ll 
for each of the detector heights 3, 10, 20, 50, and 100 feet. If the area 
does not meet the criterion, it should be subdivided into areas which do. 


C. Shielding 

In the previous sections we have considered only unshielded detectors. In 
this section we shall obtain conventional multiplicative factors to account for 
shielding. The cases considered are: 

(a) horizontal shielding between the detector and a centered area 
of contamination; 

(b) horizontal shielding between the detector and an off-center area 
of contamination; and 

(c) vertical shielding between the detector and an off-center area of 
contamination. 

1. Centered Areas and Horizontal Barriers 

The horizontal barrier shielding factor for centered areas B q is given 
by the equation 


B (X Jh-z|,A) 
o o ■ ' 


C Q <X o ,|h-t|,A) C 0 (X o . U ) 
C o (0,|h-z|,A) “ C o (0,W) 


(D-ll) 


where 


X » shield equivalent mass thickness in psf 

|h-z| » vertical separation of detector and source in feet 
A “ area in square feet 

w ■ solid angle fraction subtended by area A at 
height j h-z j 


Co <X o , ) 


overhead contribution for the given X and w 
(Reference D-3) ° 
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Figure D-ll 

Maximum Square Area To Be Represented by a Point 
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Horizontal Farrier Factors For Centered Areas 
















Since cMX^u*) is a function only of the solid angle subtended, separate 
graphs are not given for each |h-z| in Figure D-12. Rather, a nomogram is given 
in Figure D-12 for converting areas at different |h-z| values to equivalent areas 
at jh-z| * 10 feet. Essentially this is the same approach as that taken in 
Reference D-5. 

The B q (X o> |h-zj,A) values of Figure D-12 are used in the following manner 
to obtain the shielded dose contribution from a centered area: 

(a) find C d (jh-zj,A) using Figure D-18 

(b) in the Figure D-12 nomogram locate the given area value on the A 
scale (on the left) 

(c) draw a straight line from this point A''through the given value of 
{h-z | on the middle scale and continue to the A''scale on the right 
side 

(d) from the intercept of thi9 line with the A'' scale, follow a horizontal 

line to tne curve for the correct X (if X falls between two curves, 
o o 

an interpolation should be made) 

(e) read the B (X ,jh-z|,A) value at the bottom 

o o 1 1 

(£) multiply C d (|h-zj,A) by B q (X o ,|h-z|,A) to obtain C d (X q ,| h-z|,A). 

This would be the contribution from the centered area Areas also 
shielded by vertical barriers should be treated as off-center areas. 


2. Off-center Areas and Horizontal Barriers 

The horizontal barrier shielding factor for off-center areas B Q1 
by the equation 


(X,,jh-*j,r) 


C o (X^jh-zj,*(«*&■ - C o (X 0 ,|h- 8 |,*(*-&•)-) 

(0, |h-*j ,s(r4^^) - C o (0,(h-z: 


is given 


(D-12) 


where the terms ere the same as those defined for Equation D-ll, except that the 
area has been give*. by a function of the horizontal separation of aource and 
detector. 
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The values of B q (XQ,|h-z|,r) are given in Figure D-13. Again, separate 
curves are not given for each |h-z|. The C q (X Q ,|h-zj,A) values are functions 
only of the solid angle subtended by the Area A which is defined by r in 
Equation D-12. The nomogram thus can be used for converting radii at different 
|h-z| values to equivalent radii at |h-z| « 10 feet. The values of Figure 
D-13 are used in analogous manner to those of Figure D-12. That is, to obtain 
the shielding dose contribution from an area: 

(a) find C^(|h-zi,r) from Figure D-8 and multiply the value by the area 
(a square area satisfying the criteria of Figure D-ll) 

(b) in the Figure D-13 nomograms, locate the given horizontal separation 
distance r on the r scale (on the left) 

(c) draw a straight line from this point r through the given value 
of |h-z| on the middle scale and continue to the r* scale on the 
right side 

(e) read the (X Q ,|h-zj,r) value at the bottom 

(f) multiply [C, (|h-zj,r) J A by B Q1 (X Q , |h-z| ,r). (If there are no 
vertical barriers, this is the contribution from that area.) 

3. Vertical Barriers 

As a nultiplicative barrier factor for vertical barriers, we again need 
the ratio of the contribution with the barrier present to tbs contribution 
without the barrier. First, the case for which no horizontal barriers exist 
between the source and the vertical barrier is considered. The contributions 
through the walls are obtained from values in "The Equivalent Building Method" 
(Reference D-5). To obtain these contributions, we find the protection factors 
(PF) with no roof contribution (infinite root weight) and we note that the P? 
is the reciprocal of the relative contribution through the wall. Thus we 
obtain 
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Horizontal Barrier Factors for Off-Center Areas 






























(*) 


C d (X v ,A') PF (co,0.A ) 

C d (0,A*) PF (»,X v A V ) 


(D-13) 


where 


B (X ) » vertical barrier factor with no horizontal shield 

v v 

between the source and tl ? vertical barrier 

X = equivalent mass thickness, in psf, of the vertical barrier 

* * 2 
A - (2r y 

r * horizontal distance from detector to the vertical barrier t 
The PF values in Reference D-5 are given for creas A* in values ranging from 
100 to 100,000 square feet. Since these numbers include wall scattering as 
well as direct attenuation, it is conceivable tliat the barrier factor (X^) 

is also a function of the distance r . However, the ratio of equation D-13, 

* 

when obtained for the complete range of areas A given in Reference D-5, shews 
nc appreciable difference between the values obtained for tie different areas. 
That is, the variance introduced by reading the graphs of Reference D-5 was 
greater than the differences in the values for different, areas. Thus, a 
single curve of B^ (X^) is given in Figure D-14. 

If a horizontal barrier exists between the source and the vertical 
barrier, the curve derived from Reference D-3 for this configuration will be 
used. This curve is reproduced in Figure D-l' as (X^). 

If there are apertures in a vertical barrier, corrections to the mass 
thickness of the solid portion should be made before using Figure D-14. In 
the configuration with an intervening horizontal barrier, the mass thickness 
of the solid wall should be multiplied by the fraction of the wall not 
voided by apertures. Thus, 

X « (1-Ap) X (D-14) 
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Figure D-14 

Vertical Barrier Shielding Factors 
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where 


Ap * fraction of wall voided by apertures 
X * mass thickness of solid portion of wall 
In the configuration with no intervening horizontal barriers, correction 
curves for apertures are presented in Reference D-5. However, a consideration 
of their basis indicates that the correction curves of Reference D-5 deviate 
somewhat from the form which they should possess (see Chapter III). Because 

of this deviation, and since the corrections do not differ greatly for 
* 

different areas A , a single "consensus" set of correction curves will be 
used in this model. Figure D-15 presents the aperture corrections to be 
applied in this configuration with no intervening horizontal barrier. To 
use these correction curves, compute the precent of the vertical barrier 
occupied by apertures 7, Ap, and determine the mass thickness of the solid 
portion of the barrier X . For these values of % Ap and X , read from 
Figure D-15 the value of X^ to be used in Figure D-14 to obtain (X v >. 
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200 , 
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Correction to Vertical Barriers for Apertures 






































































III. EVALUATION OF EQUIVALENT BUILDING METHOD APERTURE CORRECTIONS 


4 


Shelter Design and Analysis Volume II presents the Equivalent Building 
Method of Fallout Radiation Shielding Analysis Design. The Equivalent Building 
Method is based on the approximation of reducing a complex shielding situation to 
an equivalent simple, solid-wall, single-story structure problem. For above the 
ground, the basic structure assumed in TR-20, Volume II is a solid-wall, single- 
stc y, square building. The wall height is 13 feet with the detector 3 feet above 
the floor. The aperture sill heights are 3 feet, the detector height. Basic 
structure areas of 100, 1,000, 10,000, and 100,000 square feet are used. For 
each of these areas a family of curves which relate the protection factor to wall 
mass thickness, with roof mass thickness as a parameter, is given. When a structure 
varies from the basic case, an equivalent mass thickness is found from given curves 
such that a standard structure with this wall weight has the same protection factor 
as the structure under study. In particular, corrections for apertures are made 
using Figure 10 of TR-20, Volume IX, reproduced here as Figure D-16. 

The purpose of this appendix is to examine the peculiarities exhibited by 
these curves at mass thicknesses less than 100 psf. In Figure D-16, Ap is the 
fraction of the wall taken by apertures; so let Ap be the fraction of wall above 
the detector taken by apertures, and let all apertures be above detector height. 

Only the ground contribution C g is affected by apertures in the exterior walls; 
therefore, equivalent wall mass thickness is found by equating ground contributions 

c e <ap, x e ) ■ c g (0, x;> 

Since 

(0,X ) 

a e < o 

Tx 

e 

(Figures 1 through 4 of References D-? and D-5) 
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Apperture Correct: 







c g (Ap > V < c (o, x e ) <s«o x; > x e 


(D-15) 


Using the method of the Engineering Manual (TR-20, Volume I), we have 
C g (Ap, X e ) - jtG^) (1-Ap) G a (O ) U-S u (X e )] 

+ id-Ap) <=„<«„> + = s (^) ] S u (X e )E | B c (X c ,H) 

+ *P C a <- u > * 

Given Ap < 1, and Ap <1, then by equation D-16 


C (Ap , X ) < C (Ap X ) is equivalent to 
g e — g e 


{- ^ = a (%) [l-sj - ^ G s ( % ) SE | B e + 4 G a <^> < 

{- A P G a <«U> [1-S u j - A P G s <" u > S „ E } B e + A P G a S' 


** { G .H>ii-B e *v e )- G .S> s , E » e }> 

G . <“„> t l - »„ + W * V“u> S w EB . ' 


(D-17) 


Ap ■ 0 , then by equation D-17, C (Ap , X ) < C (0,X ) if and only if 
r g r c — g e 

0 > G (u> ) 1 - B + S B - C ( u ) S EB 

— a u I e w cj a u w e 


or G (w..) 


(D-18) 


-16) 
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Now, for any Ap < Ap , satisfaction of condition D-18 implies satisfaction of 
condition D-17. Thus, with Ap >Ap, if all apertures are above the detector, then 


X" > X <=ss> C (Ap , X ) < C (0,X ) <•**'> C (Ap, X ) < C (Ap,X ) 
e - e g 1 v e y - g e' g e' - g e' 

<®> C g (0, X") < C g (0, X^) <=C> X" > X^ . 

For a given area, then X^ plotted to X g with Ap as a parameter has the form of 
Figure D-17. 



Correct Form for Aperture Corrections 

The crossover point, X @ , is given by the equality of condition D-18. In 
condition D-18, for a given eccentricity ratio, the right side is independent of 
area and the left side is independent of X^. The crossover point may then be 
found as in Figure D-18. 
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Deterninatlon of Points of Equality of Condition D-18 








IV. AN EXAMPLE ANALYSIS 


The application of the point source model is illustrated by a simple, but 
fairly realistic example. Consider the decontamination situation of Figure D-19. 

In this example the detector is located cee feet above an open area (a parking 
lot or park, say) which covers a city block. Area I includes the parking lot and 
the adjoining streets. Area II is the flat roof of a building one block long, 75 
feet wide, and 15 feet high. Area III is another parking lot separated from the 
detector by the building. 

The detailed analysis or these three areas is shown in Table D-I. First, the 
centered area is analysed; the analysis is very simple in this example since the 
eccentricity meets the criterion and there is no shielding. If this area does not 
meet the eccentricity criterion, it should b* subdivided and the off-center sub¬ 
divisions should be treated as other off-center areas,. 

Since Areas II and III are quite elongated and the areas do not meet the 
criterion of Figure D-ll, they are subdivided into areas which are nearly square 
(in this case exactly square) and which meet the area criterion in Figure D-ll. 

The data for Area II with |h-z| * 12 feet is interpolated between the 10 and 20 
foot curves. Note that the horizontal barrier factor is off-scale on Figure D-J2 
and thus is carried as an inequality. Note also that the aperture correction for 
vertical barriers is given by ■ (1 - ) X for Area II where a horizontal 

barrier exists between the vertical barrier and the source, and is given by 
Figure D-15 for Area III where there is no horizontal barrier. 

In the application of the point-source model, a tabulation sheet such as 
Table D-I will be found very useful. 
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V. THE EQUIVALENT PLANES METHOD 


The Equivalent Planes Method is based on the analyses of the previous sections. 
The keys to the method are two limiting criteria. The first gives limits on the 
length-to-width ratio (eccentricity) of an area centered above or below the detector 
such that the dose contribution is within ten percent of the approximating dose 
from a centered square of the same area. The other criterion gives limits on the 
area of an off-center square such that the dose contribution is within ten percent 
of the approximation given by the product of the area times the dose contribution 
per unit area from the center of the square. 

Shielding factors are presented for above ground detectors which are separated 
from source areas by relatively light shields. These data should be most applica¬ 
ble to decontamination studies. The dose contributions calculated are suitable 
for inclusion in formulas for decontamination analysis. 

A concise manuel, giving a step-by-step procedure for applying the Equivalent 
Planes Method, is included at the end of this volume. Several data tabulation 
sheets to be used with the manuel are also included. 
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VI. CONCLUSIONS AND RECOMMENDATIONS 


It is concluded that the point-source model can be a useful and realistic 
method of analyzing fallout radiation contributions for decontamination studies. 

The real keys to the method are the limiting criteria given in Figures D-7 and D-ll. 
The first criterion, Figure D-7, gives limits on the length-to-width ratio 
(eccentricity) of an area centered above or below the detector such that the dose 
contribution is within ten percent of the approximating dose from a centered square 
of the same area. The other criterion, Figure D-ll, gives limits on the area of 
an off-center square such that the dose contribution is within ten percent of the 
approximation given by the product of the area times the dose per unit area at the 
center of the square. These criteria are based on unshielded areas. In future 
studies it would be desirable to investigate these criteria for shielded areas. 

Shielding factors employed in this model are obtained from well-known sources 
(References D-l, D-2, D-4 and D-5). The approach taken is based on the premise that 
above ground detectors which are separated from source areas by relatively light 
shields are of primary interest in decontamination studies. If applications 
indicate that basement areas and heavier shields are of more interest, the model 
can readily be extended. Probably it will be found necessary to generate data 
for Xq values between zero and 25 psf to fill in Figure D-13 since many "curtain" 
walls and light roofs have mass thicknesses in this range. The data are not avail¬ 
able from Reference D-3, so they will have to be generated directly from Spencer's 
data, Reference D-l. 
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Appendix E 


Note: 


A FORTRAN Program for Decontamination Analysis 


The material in this Appendix was originally submitted to 
Research Memorandum RM OU-214-10.* 


USNRDL as 


C. Dillard and J. Ryan. A FORTRAN Program for Decontamination Analysis 
RM-OU 214-10. Durham, North Carolina: Research Triangle Institute, 
Operations Research and Economics Division, 31 August 1965. 
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Appendix E 

A FORTRAN Program for Decontamination Analysis 


I. INTRODUCTION 


A. General 

This appendix describes a computer program written in FORTRAIi lo compute the 
effectiveness parameters used to analyze municipal decontamination. The program 
was written in FORTRAN 64 to be used on large scale computers such as the CDC 3600. 
The program has been debugged and tested and has b e en used under OCD work unit 3233B. 

For a given level of decontamination the program is capable of determining 
both the reductions in dose-rate at specified detector locations and the re 
ductions in total dose for persons spending prescribed amounts of time at 
specified detector locations. Other parameters used in the analysis of municipal 
decontamination, such as the PF's at the detector locations and the equivalent 
PF's* associated with the activities (without decontamination), are also computed. 

B. Contents 

Section II of this Appendix is a description of the computer program and 
includes: 

1. a list of the equations which are computed 

2. the data input formats 

3. the printed output formats 

The basic equation defining the equivalent protection factor Is: 

EPF • fj/p, * ( 2 /Pj'+ 

where f. is the fraction of time spent with protection P,. 



the detailed flow charts 

a card listing of the program with a set of sample inputs 
a sample output of the program. 
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II. THE COMPUTER PROGRAM 


A. Introduction 


The CDC 3600 FORTRAN program computes the values of the following: the pro¬ 
tection factor (PF) at each location, the equivalent protection factor (EPF) for 
each activity pattern, the fraction (CF) of the total intensity prior to decontami¬ 
nation at a given detector due to a particular contaminated plane, the intensity 
reduction factor (RN) for each decontamination strategy at each detector location, 
and the activity dose reduction factor (RN^) for each combination of activity and 
deconi, .inution strategy at all detector locations. 

Calculations are made from the following input data: Contributions (C) to 
intensity at each detector location from each plane of contamination the fraction 
(f) of time that an individual spends at each detector location performing each 
particular activity, the mass reduction factor (f) for each decontamination strategy, 
and the key (S) to associate the appropriate value of F with each composite decon¬ 
tamination strategy for each plane. 

The program is designed to accept a maximum of twelve each: detector locations, 
planes of contamination, activity patterns, composite decontamination strategies; 
and twenty mass reduction factor values. 

It.j equations used to compute the above parameters are the following: 


PF 


j 



(E 1) 


EPF 


k 



CF 


ij 


c 4 

P 

i 1 


JL. 


(E-2) 


<K 3) 


E-3 


(E 4) 


RN 


jh 


"^h 



(E-5) 


where: 

j identifies the detector location; d, the number of detector locations 
j ■ 1, 2, ... , d (d meximum ■ 20) 

i identifies the plane of contamination; p, the number of planes of contamination 
i ■ 1, 2, p (p maximum * 12) 

k identifies the activity pattern; a, the number of activity patterns 
k=l, 2, a (a maximum = 20) 

l identifies the decontamination strategy; s, the number of decontamination 
strategies. £ * 1, 2, s (s maximum * 20) 

h identifies the composite decontamination strategy; c, the number of composite 
decontamination strategies, h » 1, 2.c (c maximum = 100) 

C . identifies the contribution to intensity at , detector location from ith 
plane of contamination 

f^, identifies the fraction of time spen^ at detector location j performing 
the activity denoted !>y k 

F, identifies the mass reduction factor (fraction of fallout remaining) for 
decontamination strategy ( 

3^ identifies whi^h strategy £ is used to decontaminate plane i for composite 
strategy h 

B. Data Input Format 

The first input card, identified by a 2 in card column one, specifies the 
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number of values to be entered for each variable subscript in integer form (XX) as 
9hown below: 


cc cc cc cc cc cc 

l _ 2,3 4,5 6,7 8.9 10.11 

/l XX XX XX XX XX 

e n dp a sc 


Each value for each of the input variables C., and f, . is 

ij kj 

card and entered in sequence by row (A^ A^ j*-*- A ] n > 

A* ;A -, A „, A _, ... A ). 

2,n m,l m,z m,3 m,n 


entered on a separate 
A 2,l' *2,2’ A 2,3’''' 


CC cc cc 

l _2^10_ 1 

/T xxx}ixxxx /X 

rd C . (Real Number)! *d 

1 C n 


CC 

- 

X . X X 

f^j (Real Number) 


The values for input variables (up to twenty) are entered from a single 
card. The first value, which serves also as the identification number, in card 
columns 1-4 must always be 1.00. 









C. Printed Output Format 


D = XX, P = XX, A = XX, S 

- XX, C * XX 


Values of C(J,I)-- 



Detector 


PLANE NUMBER 


Location 

1 2 

3 

p 

No. 1 

.XXXXXXXX .XXXXXXXX 

.XXXXXXXX . . . 

•XXXXXXXX 

No. 2 

.XXXXXXXX .XXXXXXXX 

.XXXXXXXX . . . 

.XXXXXXXX 

No. 3 

.XXXXXXXX .XXXXXXXX 

.XXXXXXXX . . . 

.XXXXXXXX 

No. d 

.XXXXXXXX .XXXXXXXX 

.XXXXXXXX . . . 

.XXXXXXXX 

Values of FR(J,K) — 



Detector 

ACTIVITY PATTERN NUMBER 


Location 

1 2 

3 

a 

No. 1 

X.XX X.XX 

X.XX . . . 

X.XX 

No. 2 

X.XX X.XX 

X.XX . . . 

X.XX 

No. 3 

X.XX X.XX 

X.XX . . . 

X.XX 

No. d 

X.XX X.XX 

X.XX . . . 

X.XX 

Values of S(H,I)-- 




X.XX X.XX 

X.XX . . . 

X.XX 

Values of S(H,I)-- 



Composite 

PLANE NUMBER 


Strategy 

1 2 

3 

p 

No. 1 

XX XX 

XX ... 

XX 

No. 2 

XX XX 

XX ... 

XX 

No. 3 

XX XX 

XX ... 

XX 

No. c 

XX XX 

XX ... 

XX 







PF(J) = 





Detector 

Location 





No. 1 

xxxx.xx 




No. 2 

xxxx.xx 




No. d 

xxxx.xx 




EPF(K) = 





Activity 

Pattern 





No. 1 

xxxx.xx 




No. 2 

xxxx.xx 




No. a 

xxxx.xx 




CF(J,I) = 





Detector 


PUNE NUMBER 



Location 

1 2 

3 4 

P 


No. 1 

x.xx x.xx 

X.XX X.XX 

. . . X.XX 


No. 2 

x.xx x.xx 

X.XX X.XX 

. . . x.xx 


No. 3 

x.xx x.xx 

X.XX X.XX 

. . . x.xx 


No. d 

x.xx x.xx 

X.XX X.XX 

. . . x.xx 



Check by summing 

CF(I ,J) on each J. 

Should equal L. 



X.XX X.XX 

X.XX X.XX 

X.XX 
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RN(H,J) = 

Composite. 

Strategy 

1 

2 

DETECTOR LOCATION NUMBER 

3 4 

d 

No. 1 

x.xx 

x.xx 

X.XX 

x.xx . . 

X.XX 

No. 2 

x.xx 

x.xx 

X.XX 

x.xx . . 

X.XX 

No. 3 

x.xx 

x.xx 

X.XX 

x.xx . . 

X.XX 

No. c 

x.xx 

x.xx 

X.XX 

x.xx . . 

x.xx 

RNA(H,K) = 

Composite 

Strategy 

1 

2 

ACTIVITY 

3 

PATTERN NUMBEP 

4 

a 

No. 1 

x.xx 

x.xx 

X.XX 

X.XX . . 

X.XX 

No. 2 

x.xx 

x.xx 

x.xx 

X.XX . . 

X.XX 

No. 3 

x.xx 

x.xx 

x.xx 

X.XX . . 

. „ X.XX 

No. c 

x.xx 

x.xx 

x.xx 

x.xx . . 

. x.xx 




D. 
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K-ll 




























E-14 












B-15 
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SCFI(J) *0. 













s-n 
























8 - 2 } 
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E. Card Listing of Sample Input Deck 


•SEQUENCE*339 ." 

•J03*12345*JRYAN 
•FTN.L.X 

PROGRAM JRYAN 

DIMENSION C( 10*20) «FR ( 20.20) .F (20 ) .LS( 1 0 . 10 ) * PF ( 20) .EPM2Q) ,CF<cO 

_1 * 10).RN(100*20).RNA(100.20).SCFI(20)*NUM(100) 

PRINT 1000 " . 

10CC FORMAT(50X16H GLOSSARY) 

PRINT 1001 

1OC1 FORMAJ{*OC( I .J) . DENOTES CONTRIBUTIONS TO INTENSITY AT JTH DETECTS!. 
1 LOCATION FROM ITH PLANE OF CONTAMINATION*) 

_ PRINT 100 2 __ 

1002 FORMAT(*OFr”(K *J)DENOTES THE FRACTION OF TIME*SPENT AT DETECTOR J 
2ERF0RMJ.NG THE ACTIVITY DENOTED BY K*) 

PRINT 1003 

1003 FORMAT(*0F(L) DENOTES THE MASS REDUCTION FACTOR (FRACTION CF FALLS 
3UT REMAINING) FOR DECONTAMINATION STRATEGY L* > 

PRINT _ 1004_ .. _ _ 

13~4 ‘FORMAT (*OS'(H« T)‘ DENOTES'' WHICH"STRATEGY"’L" IS USED TO OECONTAMINA ! i 
4 PLANE I FOR COMPOSITE STRATEGY H*) 

PRINT 1005 

1005 FORMAT(*0PFCJ) "NOTES THE PROTECTION FACTOR AT EACH DETECTOR LC. - 
5TI ON*) 

_ P RIN T 1C 06____ 

l0G6* FORMAT t*CEPF(io “DENOTES THE EQUIVALENT'' 'pROTECTl Cn' FACTOR FOR L A C ! - 
6ACTJ V i TY PATTERN*) 

PRINT*1007 

1007 FORMAT(*GCF(J ! _I ) DENOTES THE FRACTION OF THE TOTAL INTENSITY HF I. • 
7 TO DECONTAMINATION AT”A GIVEN DETECTOR*) 

_PRI NT 100 8........ 

Tooe FORMAT < * DUE TO a' PARTICULAR CONTAMINATED PLANE*")’ 

PRINT .1009 

1 C39 FORMAT {*ORN (H * J) DENOTES THE INTENSITY REDUCTION FACTOR FOR EAC-i. 
8DEC ONT AM I NAT I ON STRATEGY AT EACH DETECTOR LOCATION*) 

PRINT 1010 

1010 FOR MAT< » ORNA ( H.K ) DENOTES THE ACTIV ITY DOSE REDUCTION FACTOR FC- 
’*9ACh"cOM81 NAT ION OF* ACTIVITY AND DECONTAMINATION*) 

PRINT 1011 

1011 FORMAT(* STRATEGY AT EACH uETECTCR LOCATION** 

99 READ 1«ICK,ID.IP.I A,IS.1C 

1 F 3RMAT(11.512) 

- - 1F(E0F,60)41. 42___ _ 

42 IF(ICK-2>61.43,61 

43 DO 10 !•! .IP 
00 !6 J*1.ID 

REAO 2tICK.CtI,J) 

2 FORMAT * 'll »F9.8) 

_ IF( !C*-3>62.10.62 ___ 

io’continue 

DO 20^ l . I A 
DO 20 J«f. ID 

_READ 3,IC<,FR(K,J>_ _ 

5 FORMAT ( l l •>*•27 

SF ( ICK»4 >63,20.63 .... . 
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v 20 CONTI Nut 

READ 4.<F!L)«L=1.IS) 

|: 4 FORMATC20F4*2) 

f IF(F{1>-l,)64.30.64 

! 30 00 40 M=1.IC ' 

I . R£ A D 5. ICKtNUMC M ) . <LS (M. I ) ♦ I = 1 jJP)_ 

1F ( ICK—5 ) 65 « 40 .65 

40 CONTINUE_ 

PRINT 21 * ID. IP. I At isV'fc 

cl FORMAT i4H1D =.I3.3X,3HP =.I3.3X.3HA S.13.JX.2HS *.I3.3X.3HC =«l j,/1 

PRINT 22 .. .... 

_££ FORMAT ( /UPVALUES OF C ( I . J ) —*_>___ _ 

*PR1*NT 220 

FORMAT(9H0DETECT0R »53X«22HP LANE NUMBER) 

PRINT 221«ClTl = l♦IP) . 

! <£•« 1 FORMAT! 16H LOCATION .13.10111) 

DO 222 J= 1 « 1D 

*22.. PR I NT.. 1 1 t J . (Ct 1 . J ) . 1 a 1 . IP ) ; .... 

11 FORMAT (4KONO* *1 St 2X111F1 US)”* 

PRINT 23 

23 FORMAT(/*0VALUES OF FRIJ.K)—*) 

PRINT 230 

"23C FORMAT(9HCDETECTCR.42X.43HA C T 1 V I T Y PA T T E R N NUN c 

_ IE R)___.___ _ _ 

PRINT '23ft < I . 1=7 ."i A) 

23 1 FORMAT! UH LOCATION .2016) 

DO 232 J=T.ID 

232 PRINT 12. J_t (FR!K. J) »K=1 . I A) 

12 FORMAT!4h6noV«15.3X.20F6.2) 

PRINT 24___ 

24 FORMAT!/*6VALUES"OF F!L>--*/) 

PRINT 13.<fCL).Lai.IS) 

13 FORMAT!12X.20F6.2) 

PRINT 25 

c5 FORMAT!/*0VALUES OF S(H*I)—*) 

_ PRINT 250_____ 

2~0 FORMAT*(~l 0HOCOMPOS 1 TE152X~«22HP L A N E~ N U M~ 6 E R ). 

PRINT 251 «!I.1=I.’P) 

231 FORMAT!I2H STRATEGY ' .1216) 

DO 352 M« U IC .. 

252 PRINT J4,NUM!M i , !L3{M. I ) . I-' 1 . iPi 

14 FORM AT!4HONO+.15.3X.1216) __ 

do 60'j=*r.TcT~' 

SUMCIJsC. 

DO 50 1*1. IP ' 
o 0 SUMC i J * 30 MC U*Ct It _) 

6C PF ! J).« 1 t/SuMC 1J 

.. P RIN T IS_ _ _;__ 

15 FORMAT !*iPF.(J> **)’ 

■ PRINT 150 ' 

IOC r ORM AT <9H00ETECT0R. /9H LOCATION)' 

PRINT 151 »J 3.PFJU). J*1 » IQ| 

. i'l r OR«AT v*f.0N0i"tT?tF9T2» * * ■ .. 

DO 60 <= 1 « IA 


1*26 









SQ = G • 

j ....DO. 70 J=|«ip._ _ 

70 SQ=SQ+FR<K,J)/PF<J) . 

8 0 E PF (K)»U /SQ_____ 

PRINT 16 

16 FORMAT (*OEPFtK> = ») __ _ 

PRINT 160 ~ .. ’ ... .. 

160 FORMAT < 9H0ACT I VI.TY./8H PATTERN) 

PRINT 161 * (K . EPF( K)* K=1«I A > . 

161 F0RMAT(4HGN0*,I5.F9.2) 

DO 110 J=1.ID 

_ SCI J = 0 ___ __ 

DO 100 I«1* IP 
100 _SCI J*SCI J+C< i .J_)_ 

.“do” lio 1*iTip.. 

110 CF( UJl-Cl I .Jl/SCIJ __ 

.PRINT 17 ..... .. 

_ 17 FORMAT( /»0C F(J.I)= » )_ 

PRINT 220 .~ *’ ' .——-.- 

PR I NT 231.(1 * I_s 1_» I P_) _ 

DO 120‘J=1 . ID . 

120 PR I NI_ 1.2 • J » (CF ( I «_J) « I * 1 .IP) 

PRiNT 9 

_9 FORMAT(*0 CHECK BY SUMMING CF(J«1 ) ON EACH J. SHOULD EQUAL 1•*/ 

DO 130 J=ltID 

SCF I(J) =0 _ _ 

DO 130 1*1.IP 

1 30_SCFIJ_JJ = SC FIJ JJ+CFiI «J_)_ 

PRINT 13.C SCFI{J)*3=1.10) 

OQ 140 M= 1 . IC ___ 

DO 140 J=1 . ID 

_RNU).M)=0.____ 

"66 140 fs f.IP 

L = LS(M. I ) _ 

14C RnTJ«M)‘=RN(‘j7mY+CF( I . J)*F(L> 

_ PRINT 18 __ 

18 FORMAT(/#QRN(H.J) =*) " . 

PRINT 130 

ISO" FORMAT( lOHOCONiPOS I T E » 38X .45HD g > E C T 0 ft LOCATION N 

1MJ3_E R)_ ___ 

PRINT 254« (J » J=1,ID) 

_ DO 253 Msl.IC __ 

'253 PRINT 12.NUMIM),(RN<J.M),J*1«ID) 

_._D0_ 31 0 K* 1 . IA __ __ 

" DO” 310 Mil* 1C 

SUSO»0. _ _ __ 

susp =Td.' 

_ DO 300 J 3 1 . ID ___ _ _ _ 

SUSQ*SUSQ+FR1K.JJ/PF( J ) 

300 SUSPsSUSP+FR( K .J)/ PF(J) »RN (J.M) .. 

316"RNA(K7m ) *SUSP/SUSO 

PRI NT 1 9_____ 

19 FORMAT (>*CRNA(H.K)* • { 

PR INT 190_ 
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190 "FORMAT ( | vHOCcViPO?. ITF«36x’«'a3HA C T I V ! T Y P A T T i. P 

13 ER) ._ 

PR I NT 254« (K «K= 1 * IA ji. 

DO 191 M=1«IC 

1 .< 1 PR I NT 12. MUM (M) * ( RNA (K . M) . K = 1 . I A > 

_PRJ.NT__9S _ __ 

98 FORMA tTihI ) 

GO TO 99 

61 PRINT 71 

71 FORMAT!* 0 INCORRECT DAJA_CARp—-NOT VALUES OF SUBSCRIPTS.*) 
" GO TO'41 

62 PRINT 72 _ _ 

72 "FORMAT! *6*1 INCORRECT DATA CARD—-NOT VALUE OF CU.J).*) 

GO TO 41 

c3 PRINT 73 

73 FORMAT(*0INCORRECT DATA CARD—NOT VALUE OF FRIK.J).*) 

GO TO 41 

-i4 PR IN T _7*.____ 

74 FORMAT f*0 FnCORRECT DATA “CARD—NOT* VALUES OF* F <L >"• *) 

GO TO 41 

o5 PRINT 75 

75 FORMAT(*0tNC0RRECT DATA CARD—NOT VALUES OF S<H.Ii.*) 

GO TO 41 

. 5..F0RMAT.CJ J_2I g L.___ 

<-54 FORMAT (1 1H STRATEGY .1216*)" 

41 CONTINUE 
END 

SCOPE 

•LOAD 

♦ RUN«3.2CCC . 1______ 

20404050606 

5.299 

2.336 

.024 

3.002 

3.031. *. . . ~ 

3.2 09 
3.005 
; • _ 29 
U206 
• 1 82 


-..071 

,,j77 

j.085 

-.002 

4-0.06 

-0.07" 


- ..20 
-C .67 



i-2 r. 









• f- 7 
0 »jr.fi 
0.02 
J . 0 3 
..67 

0.05.. 

u. 15 
w.67 
l.. 19 
o.C8 

o.06 

0.67 

.000.020.030.090.100•13 
:'0 1C 1050102 
:’?r303C401 
0301010401 
.0401030601 
'■ ; •'.0 4020 101 
w0601010505 
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F. Sample Output 


l"2?473 JM=A3 1440? 

1 00213 IU». 0:347 

137411 TSm. 

»»03»*h fciitms. 

•tour 

LABELED C3«40M 

*DMRfcRr.D :0«404 
NON* 


1 63050 ioS. 00423 
1 5770C ALLOC, 00?i3 
1 37330 flea|7IOC 00033 


1 01040 USOEflRn* 0020? 
1 37313 RSTA4 00243 
1 33141 ION, 02213 


1 0l3A2 flOOEXTAV 00004 
1 37403 STM, 000*0 


LOTOV 031VTS 
4 000»4 SiXlMf 
i 01302 aofloi: t, 
1 02301 lOSIlSl. 
1 07242 OlflCMAlN 
1 01040 }03k»40H 

1 37001 4$CHfcC4 
1 37700 4ETJ4V, 

1 33103 IflM, 

1 37213 .0E7C37, 


3334 JH*»* 

7305 OODITEflf 
2050 100. 

2173 flOOMlST, 
1733 fliOEOSET 
7010 RSClOSE 
0072 OUST. 
3141 OCL'Bur, 


1 01300 OOOEXTAV 
1 57415 TSm. 

1 07373 100. 

1 02304 axnOJBL, 
i 01303 exit 

1 37577 XOIStltX 
1 00101 IRETURX, 
1 37301 oooinoc 


-A 


02307 TMFN0, 

mil stm, 

1 02070 IOS, 

['tint TTn»; - 

1 373l3 AST OR 
1 0011* ALLCe, 

1 00172 ALLOC IM, 

I »oW TlB.— 


fcOECUT103 STAMIEu AT 1417 >02 


BOO 













C<I,V» DENOTE* CUNTRISuTIONS TO INTENSITY *T JIM OtTECTOR L°CAl|0.. : WON |TM PLANE OP CONT AH 1 NAT I ON 

fRI«,J)OENOTES the ENACTION or TINE spent at DETECI ON J PERFORMING the activity denoieo by K 

Fltl DENOTES The MASS REDUCTION FACTOR (FRACTION OF FALLOUT REMAINING) FOR DECUNlAMI NAT I ON STRATEGY L 

SIM,|> DENOTES nhICM STRATEGY L IS USED TO DECONTAMINATE PLANE I FOR COMPOSITE STRATEGY N 

PFlJI DENOTES THE PROTECTION FACTOR AT FaCN UETFCTOR LOCATOR 

EPFIR) OENOTES TnE EGJIVALENT PROTECTION FaCIOR FOR EACH ACTIVITY PATTERN 

CFIj.I) DENOTES the FRACTION OF THfe TOTAL INTENSITY PRIOR TO DECONT ahInaTION AT A GIYEN DETECTOR 
DUE TO A PARTICULAR CONTAMINATED PLANE ~ - 

rnim.jf dsndtes ire intensity reduction factor ruR each decontamination syrateqt at each oetector location 
Rnaim.ri denotes t me activity DOSE reduction r AfI or for each cumrinaTion if activity and OfCONT AHINaTTON 

STRATEGY At EACH OtTECIOR LOCATION 
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0 a 4 * • 4 A « 9 S « 6 C a '• 

values or cd,j>~ 

dbtecton plane n cnri nr 

LOCATION 1234 

NO, 1 0,29900000 0,10100000 0,92900090 0.07100000 

NO, 2 0,30600000 0,03100000 0,20600000 0,37/00000 

NO, 3 0.02400000 0.20900000 0.10200000 0.08*00000 

NO, « 0.00200000 0.00900000 0.00100300 0.00200000 

VALUES Of FNCJ.iO — 

obtecton activity p AXTTrmnrmr 

LOCATION 12349 

NO, 1 0,06 0.0> 0,26 0.13 0.19 

NO, 2 0.0/ 0.17 0,02 0.09 0.00 

NO, S 0,20 0.09 0,03 0.19 0.06 

NO, 4 0.6/ 0.67 0,67 0.67 0.67 

VALUES OF Ml)- 

1.00 0.02 0,03 0,09 0.10 0.13 

values or sin,i 

COMPOS!TE PLANE N U M N E N 

STNaTEQV 1234 

NO, 1 19 12 

NO, 2 3 3 4 1 

NO, 1 114 1 

NO, 4 13 6 1 

NO, 9 4 2 1 1 

NO. 4 119 9 
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-OlTICTOr 

location 

NO, 1 1,00 

NO, 2 1.00 

NO, S 2,00 

NO, A 100,00 

ioriKi • 

ACTIVITY 

httim 

NO, 1 4.22 

NO, t 0.43 

NO, S 3.11 

NO, 4 3,02 

NO, 9 3.26 

CHJ,|)« 

asTecro* plane 

LOCATION 1234 

NO, 1 0,30 0.10 0,93 0.07 

NO, I 0,39 0.03 0,21 0.30 

NO. 3 0.09 0.42 0,36 0.17 

NO, 4 0.20 0,90 0,10 0.20 

cmCck av sunning cptj.n on each j. should eoual i, 

1.00 1.00 1,00 1.00 




nn(m,j> • 

COMPOSITE 

$»»atht 
NO, 1 

NO, 2 

NO, 3 

NO, A 

NO, 9 

NO, 6 


1 2 
0.04 0.60 

0.13 0.41 

0.92 0.01 

0.44 0.70 

0.63 0.A2 

0.4* S.«4 


3 4 

0,46 0.39 

0,22 0.23 

C,67 0.91 

0,20 0.43 

0,?$ 0.33 

9.N2 0.73 


deticto* location n u m aga 


E-33 












COMPOSITE 

STRaTEOV 

NO, i 

NO, t 

NO, I 

KJ , 4 

NO, f 

NO, 6 


1 2 3 

0,99 0,63 0,00 

0,25 0.33 0,15 

0.66 0.72 0,55 

0,46 0,62 0,46 

0,56 0.60 1,62 

0,50 0.40 0,47 
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The Nature and Scope of Compand and Control System Elements 
Required for Conducting Effective Decontamination 
Analysis 


Note: The material in this Appendix was originally submitted 
to USNPJDL as Research Memorandum RM-OU-214-9+. 



Operatlens Research and Economics Division, 30 July 196S. 
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Appendix F 


The Nature and Scope of Comaand and Control System Elements Required 

for Conducting Effective Decontamination in Municipalities 

I. INTRODUCTION 


A. Purpose 

The purpose of this study is uo determine the nature and scope of the command 
and control system elements which are required to effect practical decontamination. 
Emphasis is placed on decontamination within municipalities. Such a study is 
necessary in order to answer a number of questions basic to decontamination analy¬ 
sis. The five questions on which this study focuses are: 

1. What are the preattack and postrttack data required to effect decontami¬ 
nation operations? 

2. What are the essential components of the information system needed to ef¬ 
fect decontamination operations? 

3. How should trained and untrained personnel and decontamination equipment 
be prepositioned, organized, and controlled? 

4. How can a decontamination system in a municipality be evaluated? 

5. How can a decontamination system in a municipality be most effectively 
modeled to provide a ready vehicle for system analysis? 

Of course, all of these questions are asked for various levels of decontamina¬ 
tion capabilities, requirements, and attack environments. As implied by the ques¬ 
tions above, a subsidiary purpose of this study is to develop a procedure for ana¬ 
lyzing a decontamination system in a municipality. 

B. Simms ry 

This appendix first describes elements and purposes of command and control 
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systems in general. The various components of a decontamination system in a municipality 
are then identified and embedded in the general command and control system framework. It 
is shown that a command and control system for decontamination operations must pro¬ 
vide both for decisions on •*.,<;ther or not to undertake a mission and for manpower 
and decontamination resource commitment and allocation decisions. These decision 
functions require an elaborate information subsystem consisting of organized data 
files containing prestored (preattack) data and postattack Assessments (including 
system feedback). 

The detailed characteristics of the individual components in each of the essen¬ 
tial subsystems of a decontamination command and control system are studied. The 
interrelationships among the individual subsystem components are identified and dis¬ 
played. 

In order to determine the command and control system elements required to ac¬ 
complish practical decontamination missions in a municipality, the environment and 
the system goals are reviewed and analyzed. Basic system evaluation criteria are 
also discussed and the essential decontamination system evaluation criteria are 
identified. 

Recomnendations and guides leading to the design of a basic conmand and control 
system for municipal decontamination are indicated. 




II. GENERAL CHARACTERISTICS CF COMMAND AND CONTROL SYSTEMS 


A. Introduction 

Most persons who are concerned with business, government, or military operations 
are also acquainted with the basic characteristics of command and control systems. 
However, a review of the elements and purposes of command and control systems is 
useful in providing a framework into which the components of a decontamination sys¬ 
tem can be embedded. 

Further, since the technology of command and control systems is advancing rap¬ 
idly, che terms which define the basic components of such systems are changing day 
by day. Thus, it is necessary to define precisely the terms and expressions used in 
this appendix to define the system components and system goals. 

The succeeding parts of this explanatory section are: 

Section II-B The Purpose of Command and Control Systems 
Section II-C The Basic Components of Command and Control Systems 
Section II-D Some Examples of Command and Control Systems 

B. The Purpose of Command and Control Systems 

The purpose of any command and control system is to: (1) determine for a given 
operation or set of operations, what must be done, who (or what) does it, and when, 
where, and how they do it; and (2) direct the activities which accomplish it. 

Obviously from the above, it is seen that a coesnand and control system must have 
the capacity both to initiate and to monitor activities. As operations or sets of 
operations become more complex, systematized coaoand and control mechanisms become 
more necessary. 

C. The Basic Components of Compand and Control Systems 

Most command and control systems are composed of a set of basic elements called 
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subsystems . The expressions used to identify these subsystems as well as their de¬ 
fining characteristics are presented here along the lines indicated in Reference F-l. 
The following subsystems are common to all command and control systems: 

1. Sensor subsystem; 

2. Effector subsystem; 

3. External comnunication subsystem; 

4. Internal communication subsystem; 

5. Information subsystem; 

6. Decision subsystem; 

7. On-line control subsystem; and 

8. Off-line control subsystem. 

The command and control center is composed only of the internal communication 
subsystem, the information subsystem, the decision subsystem, and the on-line control 
subsystem. Figure F-l shows the functional relationships among the basic subsystems 
of a command and control system. 

1. Sensor Sut„y3tem 

The besi'' purpose of the sensor subsystem is to acquire information inputs. 
This information input material is classified into three categories: 

a. observed information inputs; 

b. derived information inputs; and 

c. command information inputs. 

Unfortunately, sensor subsystems also acq>. re material which is not information 
input at all. This non-information input or "noise" represents irrelevant or 
erroneous material. 

Observed information inputs can be such things as visual observations, 
verbal communications, radar scope readings, radiological dosimeter readings, 
or on-ait''; damage computations. Derived information inputs are data which re¬ 
sult from deductive processes, it may be the product of either experimental 






Figure F-l Functional Diagram of Coonand and Control Syatem 
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or theoretical research, for example. This- is usually prestored data and thus 
the external sensors which acquire such information inputs may be inactive dur¬ 
ing most of the time that the system is operating. Command infonnation inputs 
are in the form of directives which enter the system from external sources. 
Decisions made outside of the command and control system but affecting the 
operation of the system are a form of command information input. Figure F-2 
shows the structure of information Inputs which arc acquired by the sensor sub¬ 
system. 

2. Effector Subsystem 

The purpose of the effector subsystem is tc» carry out: the directives of 
the decision subsystem. This may consist of simply disseminating information 
or at the other extreme carrying out the over-all system missions. The 
effectors are thus controlled by the output conmands and decisions of the 
command and control center. The effectors' task is to apply the energy com¬ 
mitted for some external action. An example of an effector would be an offen¬ 
sive aircraft directed to attack a target as prescribed in the command and 
control center. 

3. External Communications Subsystem 

The external coianunications subsystem is the network of paths which 
carries material (information inputs, infonnation, etc.) to and from the func¬ 
tional blocks tied to thfs subsystem in Figure F-l. These paths can consist 
of radio coosnunications, telephone links, verbal infonnation transmission, etc. 
Often an external sensor or effector can also be a part of the external commu¬ 
nications subsystem as well as the sensor or effector subsystem. The organi¬ 
sation and other systems variables involved in this subsystem depend largely 
on the characteristics o£ the particular coesnand and control system of which 
it is a part. 
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4. 


Internal Comcunicatlons Subsystem 


The internal communications subsystem is similar to the external communi¬ 
cations subsystem except that it is entirely contained in the command and con¬ 
trol center. Figure F-l shows its functional relationship to the system as a 
whole. 

5. Information Subsystem 

The information subsystem is the part of the command and control center 
which handles, stores, transform a, or displays the Information which is trans¬ 
mitted to the command and control center. It consists of all of the people, 
machines, reports, and files within the conmand and control center which are 
employed to perform this function. 

6. Decision Subsystem 

The decision subsystem (illustrated in Figure F-3) might alao be called the 
command subsystem. It is the component of the command and control system which 
transforms information into decisions. These decisions are usually translated 
into commands to perform an activity, or commit resources to perform prescribed 
actions during a prescribed time interval. Figure F-3 illustrates how the de¬ 
cision subsystem transforms information into decisions and comnands. 

7. On-line Control Subsystem 

Ideally, the on-line control subsystem (illustrated in Figure F-4) controls 
all information transmission and monitors all actions performed by the conmand 
and control system. In practice, the on-line control mechanism only attempts 
to perform these tasks In most real situations. Many activities are performed 
by the effectors by their own volition because there is not enough time for the 
decision subsystem to outline a course of action. For example, a Polaris sub¬ 
marine under attack presumably would either attack Its adversary or take eva¬ 
sive action immediately without communicating through a higher echelon decision 
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Figure f-3 The Declelon Sub*y* r «« •• • Tr»nafort* cion 
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subsystem. 

A very important aspect of such on-line control is the regulation of sys¬ 
tem accuracy and/or response time. This means that the on-line control sub¬ 
system attempts to control delays and errors in the transmission and trans¬ 
formation of information in the conmand and control system. This role of the 
on-line control subsystem makes it the most important element in the system. 
Unfortunately, it is also the most difficult to identify precisely in most 
physical system. In the above Polaris example, the effector subsystem per¬ 
formed a function which, for at least one situation, ioentified the effector 
as a part of the on-line control systo 1 . Usually those machines, people, etc., 
which perform on-line monitoring and control functions from within the command 
and control center are considered to comprise the on-line control subsystem. 
Figure F-4 shows schematically how the on-line control subsystem performs in 
the system environment. 

8. Off-line Control Subsystem 

The off-line control subsystem monitors and controls processes which affect 
the system but are not performed on-line. That is, they are either performed 
before the system is operational or are performed outside of the control of the 
on-line control subsystem. One important function of the off-line control 
subsystem is to collect and process statistics and research for future augmen¬ 
tation of the static variables of the system. On the basis of the observed 
statistics and research, the off-line control subsystem makes recommendations 
to the on-line control subsystem for possible future control of the dynamic 
variables. Sasic research performed long in advance of an operational system 
must be controlled, it is the purpose of the off-line control subsystem tc 
control this research. Tims, most often, scans segments of the off-line con¬ 
trol subsystem ara the first operational parts of a command and control sys¬ 
tem. 
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Figure F-4 Functional Characteristic* of On*line Control Subsystem 
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D. Some Examples of Command and Control Systems 


Some of the Air Force's "L" Systems are excellent examples of military command 
and control systems. A few are summarized below (See Reference F-2 for these and 
additional examples): 

412L - Air Weapons Control System 

An overseas theater tactical air weapons control and warning system. 
416L - SAGE Air Defense System 

A semi-automatic area air weapons control and warning system including 
the Back Up Interceptor Control (BUIC) for detecting, identifying, 
cracking, and providing interceptor weapon direction against air- 
breathing threats to the United States and Canada. 

474L - Ballistic Missile Early Warning System 

A system to provide early warning of a mass ballistic missile attack on 
the North American continent from the North. 

4811. - Postattack Conanand and Control System 

A system to enable the Conmander-in-Chief, Strategic Air Command, to 
control his forces in the event that a nuclear attack destroys or 
seriously degrades his normal facilities. 

Many of these as well as other of the *'L" Systems are subsystems of higher 

echelon command and control systems. 474L (BMEWS), described above, is also a 

sensor subsystem of higher echelon systems. 
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III. A COMMAND AND CONTROL SYSTEM FOR MUNICIPAL DECONTAMINATION 


A. General 

The previous section defined the characteristics of subsystem interrelationships 
comnon to all command and control systems. It is the purpose of this section to 
identify the basic system elements required by municipal decontamination and to 
embed these elements into the general conmand and control systems framework of the 
previous section. In order to accomplish this, the following steps are taken in 
this section: 

1. Analyze the environment, within which municipal decontamination is expected 
to operate. 

2. Identify the goals of municipal decontamination. 

3. Prescribe the decontamination system performance requirements. 

4. Identify and relate the elements of a decontamination system. 

5. Relate the elements identified in 4 above to the general conuand and 
control system framework defined in Section II. 

B. Environmental Aspects of Municipal Decontamination 

Municipal decontamination by itself is only a part of the broadsr system of 
postattack radiological defense. The time during which a radiological defense sys¬ 
tem operates can be logically divided into three time phases: the emergency phase, 
the operational recovery phase, and the final recovery phase (Reference F-3). The 
underlying reason for dividing radiological defense into three time phases is the 
change due to fallout decay. 

Figure -5 shows that the emergency phase begins with a preattack warning--an 
event highly probable during times of international crisis--if it occurs, or attack 
warning otherwise. It may last for several days, weeks, or month# depending on the 
warning obtained (combined with the population discipline attained for such an event) 
and the attack parameters, Tha decision-makers during this phase utilize preattack 
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planning, before and after warning, the extent and effectiveness of which is de¬ 
termined by the amount of preattack warning given. This is a function of intelli¬ 
gence indicators. 

The objective of radiological countermeasures during this phase is survival, for 
it is during this period that the fallout arrives, accumulates, attains some maximum 
intensity, and then begins to decay. It has been shown that the optimal or primary 
radiological countermeasure during this period is enshelterment. The central re¬ 
quirement of the system during this period is to provide a sufficient number of 
adequate shelters so located as to minimize casualties. The potential effectiveness 
of the system during this period is determined by the degree of protection afforded, 
space availability and accessibility, and warning time. 

The operational recovery phase follows the emergency phase and encompasses the 
operations prerequisite to: 

1. Sustaining life in a hostile environment; and 

2> The recovery of essential postattack activities and/or facilities. 

From Figure F-5 it is seen that the transition to the operational recovery phase 
takes place when limited egress from shelters is feasible. 

The final recovery phase begins when radiation intensity decays to an insignif¬ 
icant level and primary considerations can be focused on functional restoration of 
the target area as nearly as possible to its preattack condition. 

Municipal decontamination will be activated primarily during the operational 
recovery phase . The following environmental conditions will then prevail: 

1. Fallout deposition, generally, will be complete. 

2. Fallout radiation will still be a hazard in some areas. Moat of the 
civilian population will be sheltered during much of the time that de¬ 
contamination is being performed. 

3. There may be damage In the area from direct weapons effects (blasts, 
thermal, EMP, etc.). 
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4. The entire Civil Defence System will be very active , i.e., communications 
lines, logistics supplies and equipment, etc. will be very busy . 

Xhi3 last item is perhaps the most important condition affecting the environment 
of municipal decontamination systems . It is a factor over which the decision-maker 
has some control. As is pointed out in Reference F-3, the decision-maker has little 
or no control over most environmental factors. 

Reference F-3 also points out that the environment is an important factor lim¬ 
iting the range of decisions which can be made in the postattack period. Because 
of the wide variation of damage, yielding a wide variety of physical environments 
(in which municipal decontamination might operate), the decision subsystem of any 
proposed decontamination system must be extremely flexible and capable of making 
very complex decisions. Further, it should be capable of making decisions based on 
varying amounts of reliable (and unreliable) information. This will subsequently 
be discussed in some detail. 

C. Goals of Municipal Decontamination 

The primary goals of decontamination within a municipality (in a very broad 
sense) are: 

1. Continued survival within the municipality, 

2. Accelerated municipal recovery, and 

3. Assistance in the recovery of a region (not necessarily containing the 
municipality '.-here the decontamination system is operating). 

Since the three goals listed above are also primary goals of the overall radio¬ 
logical defense system, any decontamination system must be closely coordinated with 
other radiological defense systems. This implies the existence of an overall radio¬ 
logical defense command and control system. Although it is not the purpose of this 
report to elaborate on the overall system, it should be recognized that much of the 
cociaand Information input to the decontamination command and control system will be 
output from such a higher echelon decision subsystem. 
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D. Municipal Decontamination System Performance Criteria 

In light of the existence of a higher echelon command and control system, the 
question might be asked: "Why a command and control system for decontamination?" 
There are two basic reasons for proposing a separate system to command and control 
decontamination operations in a municipality: 

1. The specific decisions involved with effecti ig decontamination operations 
in a municipality are in themselves complex enough to warrant a "decontami¬ 
nation-level" decision subsystem, 

2. The information requirements for effecting useful decontamination decisions 
are in themselves very involved and require a separate "decontamination- 
level" information subsystem. This is not to say that a decontamination 
information subsystem would not borrow components from, or overlap with, 
parts of other radiological defense information subsystems, e.g., RADEF. 

Later in this appendix, the information input and decision structure requirements 
for effecting decontamination will be identified and will substantiate the above 
reasoning for a separate decontamination command and control system. 

The basic criteria upon which to judge a decontamination system are: 

1. To what extent does the system increase the effectiveness of decontamination 
in achieving the goals and requirements prescribed by higher echelon decision 
subsystems? 

2. How effectively does the system reduce on-line decontamination costs (rn 
manpower, crew doses, expenditures of fuel, water, etc.) for specified 
levels of performance? 

3. How effectively doss the system reduce off-line decontamination costs (re¬ 
search funds required, training of personnel, procurement of equipment, etc.) 
for specified levels of performance. 

This last criterion is separated from the second, in that it can be treated as 
a static goal of .he system (l.e., one which is evaluated by off-line simulation 
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and/or analysis of on-line system performance). The first two are dynamic goals 
(i.e., those evaluated by on-line performance and which can be controlled on-line). 

E. Elements of a Decontamination System 

The purpose of this section is to identify the functional decontamination de¬ 
cision elements as well as to determine the data base and information flow required 
by these elements. These data and information flow may be internal to the decon¬ 
tamination system (between subsystems) or external to the system (to and from higher 
echelon or para'/lei systems). The higher echelon command and control system with 
regard to decontamination is the Postattack Civil Defense Emergency Operating Sys¬ 
tem. Parallel systems would include systems to control such functions as rescue, 
law and order, engineering, welfare, firefighting, medical, etc. Certain other 
subsystems of the Emergency Operating System would provide sensor information inputs 
for the decontamination system. These include RADEF, damage assessments, NUDETS, 
etc. Figure F-6 illustrates schematically the information flow between the various 
comnand and control systems related to decontamination. 

The subsystems of the decontamination system must provide the basis for command¬ 
ing and controlling all of the activities and decisions which comprise decontamina¬ 
tion. Thus, a first step towards identifying the elements of a decontamination 
command and control system is to list all of the important activities and decisions 
which are required to perform effective decontamination. In order to determine 
whether or not to decontaminate a given facility or activity and thence to schedule 
the operation, the following steps must be taken: 

1. Determine that the facility or activity is essential to sustaining life or 
accelerating recovery. This may involve no more than being told by the 
higher echelon decision subsystem that a given facility or activity is 
essentiai. 

2. Determine that the given facility or activity is presently denied bv fall¬ 
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3. Determine that dacontamination can reduce the fallout denial time of this 
facility or activity to or below the minimum desired. 

4. Determine that equipment and supplies (water, etc.) are available to decon¬ 
taminate around this facility or activity. 

5. Determine that sufficient manpower (taking dosa histories into account) is 
available to decontaminate this facility or activity. 

6. Deternine that the facility or activity can be decontaminated to a specified 
level of effectiveness without overexposing the decontamination crews. 

7. Determine that direct weapons effects or congestion will not impede the de¬ 
contamination equipment in reaching the facility or activity. 

8. Decide to decontaminate the given facility or activity. 

9. Commit resources and schedule the decontamination of the given facility or 
activity. 

10. Initiate and conduct the decontamination operation. 

11. Monitor and control the decontamination operation. 

12. Determine whether additional resources are required and available to complete 
the decontamination of the facility or activity. 

13. Conmit additional resources to decontaminate the facility or activity if 
needed. 

14. Determine that the decontamination operation is completed. 

Each of these activities triggers a response within the system. Thus, none of 
these activities can be omitted in the analysis without reflecting a gap in the re¬ 
quired data base and organisational structure. Each activity demands a supporting 
data base, some resources, and some organisational (cooaand and control) structure. 
This organisational structure should provide for both information channels and de¬ 
cision structures. 

Figure F-7 show* the supporting structure and data bate for conducting the de¬ 
contamination of a given facility or activity. Each of the important activities 
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involved with such an operation is delineated and the data base for each activity 
specified. Clearly the required information channels comprise a complex network 
and imply the necessity for a well organized information subsystem. 

Close examination of Figure F-7 reveals the necessity for other subsystems 
similar to those described in Section II of this report. For instance, an on-lins 
control system will be required to monitor and control the decontamination operations 
and also to direct and insure that information is forwarded to the proper functional 
elements of the decontamination system. The parallelism between required decontami¬ 
nation system dements and classical command and control subsystems is considered 
next. 

F• Decontamination System Elements 

The purpose here is to identify the command and control system elements re¬ 
quired to conduct effective municipal decontamination. The interrelationships be¬ 
tween the functional elements are discussed and displayed. 

In Section II of this report, each of the necessary subsystems of any command 
and control system are identified and defined. The basic supporting data base as 
well as the activities and decisions necessary to conducting postnttaci- decontamina¬ 
tion are illustrated in Figure F-7. It is now necessary to relate these necessary 
decontamination system activities to the command and control framework described in 
Section II. Table F-I lists some of the basic dr .contamination functions performed 
by each of the command and control subsystems. Figure F-8 illustrates the functional 
relationships among these decontamination activities. It is seen that some of these 
functions are also parallel systems to decontamination. Radiological monitoring, 
for instance, would be required for almost all of the parallel systems shown in 
Figure F-b. 

What is possible to note from Table F-I and Figures F 7 and F-8, however, is that 
the s ystem required to conduct effective dr ontamination doea have all of the character¬ 
istics of a full scale command and control cystea. All of the usual aub-syatema of 
command and control -y&t<m$ are needed in some fashion to perform municipal decon¬ 
tamination. 
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TABLE F-I 


Decontamination Subsystem Functions 


Sensor Subsystem 


1. Radiological monitoring. 

2. Field Surveys, etc. 


Effector Subsystem 

1. Move decontamination equipment to sites. 

2. Perform decontamination. 


External Communications Subsystem 

1. Transmit data and information to and from radiological monitors, field 
inspection teams, etc., and decontamination crews to the Emergency 
Operating Center (EOC). 

Internal Communications Subsystem 

1. Transmit data and information within the EOC. 

Information Subsystem 

1. Store and catalog decontamination planning guides and other related 
(prestored) materials. 

2. Display environmental information (RADEF, etc.) within the EOC. 

3. Process data inputs and information coming in from the external sensors 
into a usable form for persons with information needs. 

4. Provide information to decision subsystem on request. 

Decision Subsystem 

1. Decide when, what, and how to decontaminate. 

2. Decide when to commit resources and which resources to commit. 

3. Decide whether additional information is required to make a decontami¬ 
nation decision. 

On-line Control Subsystem 

1. Monitor and control decontamination. 

2. Commit resources to decontamination tasks. 

3. Control information flow, assign tasks to radiological monitors and field 
inspection teams. This included regulating response times and information 
accuracy. 

Off- l ine Control Subsyatem 

1. Monitor and control decontamination research. 

2. Preposition and inventory decontamination personnel, equipment, and 
supplies. 

3. Organise decontamination subsystems within the EOC, 
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IV. AN INFORMATION SYSTEM (SUBSYSTEM) FOR DECONTAMINATION 


A. Introduction 

Figure F-7 in the last section depicts a very complex information-flow network 
required in the decision whether or not to conduct a particular decontamination 
operation. Certainly some control over the collecting, scoring, and disseminating 
of all of these data is necessary. It is the purpose of this section to indicate 
the basis upon which an information subsystem for decontamination could be designed. 
The steps taken towards this end are along the lines suggested in Chapter 3 of 
Reference F-3: 

1. The information input requirements are discussed. 

2. The information format Including the criteria upon which to judge the use¬ 
fulness of the information are discussed. 

3. The means for obtaining, storing and disseminating the required information 
are identified. 

The last two steps required in designing an information system (identifying 
alternative systems and choosing a best one) are not taken up in this report. 

B. Information Input Requirements for Conducting Decontamination 

Most of the information inputs shown in Figure F-7 are not required for every 
decontamination decision. For instance if no blast damage has occurred in the 
city where decontamination is being performed, specific blast-damage information 
inputs would be unnecessary. It a higher echelon decisions subsystem commands that a 
facility be decontaminated, the information inputs required to determine that the 
job is necessary are no longer relevant to the decision. Thus, there are different 
information input requirements depending on the environment and the kind of decon¬ 
tamination decision being made. The environment includes not only the attack damage 
but the overall decontamination resource situation and the population dispersal. 

(This last item is a function of the attack warning time of day of attack, pre- 
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conditioning and training, etc.) Clearly the problem of identifying all of the pre¬ 
attack and postattack data requirements for effecting municipal decontamination for 
all situations would be a mammoth undertaking, Whether or not such an effort would 
be fruitful is in itself a difficult question. This report will simply list some of 
the basic requirements which appear to be a necessary input to most decontamination 
decisions. They are: 

1. Fallout information - This includes both amounts of fallout (mass) and the 
dose-rates where the fallout is deposited. 

2. The location of essential facilities and the urgency of their recovery - 
This information may be in the form of a command from a higher echelon 
decision subsystem or prestored information in the EOC. 

3. Decontamination resource inventories and schedules - This information 
should be kept up to date in the EOC. 

4. Planning guides for conducting decontamination - This information is pre¬ 
stored in the EOC. 

These basic information requirements are presented in a general manner only; 
any additional refinement would require a complete situation analysis including a 
description of the environment and the kind of decontamination decision being made. 

C. Information Utility Criteria 

As stated in Reference F-3, the basic criteria upon which to judge information 

are: 

1. Timeliness; 

2. Pertinence; 

3. Reliability; and 

4. Comprehensibility 
1. Timeliness 

The value of radiation information particularly is very dependent on when 
it is received. This is because of the decay rate of the fallout intensity. 
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Radiation reports during the build-up phase and shortly thereafter are not 
essential for future decontamination. This Information, in fact, is best kept 
out of the decontamination information subsystem. It simply clutters the system. 
Radiation information, to be useful, must be reasonably up to date and reflect 
the amount of fallout or intensity. As the decay progresses it becomes less and 
less important that the information be processed quickly. Of course late re¬ 
ports might delay a decision to decontaminate an essential facility. 

Prestored data should be organized so that the decision-maker can obtain 
information from his files as he needs it . 

2. Pertinence 

The information should be organized and disseminated in a manner to minimize 
the amount of unnecessary information being pro ided to the decontamination de¬ 
cision-maker. Prestored information (decontamination planning guides, etc.) 
should be concise and to the point. As is pointed out in Reference F-3, "Extra¬ 
neous material not only wastes the time of those who must use it, it may also 
lead to incorrect decisions by obscuring the most important considerations." 

3. Reliability 

Information should be accurate. Reference F-3 points out that accuracy is 
influenced by the quality of the information inputs as well as the ability of 
the system to transmit and transform data without making errors. 

Unfortunately, this criterion is particularly difficult to meet for decon¬ 
tamination decisions. Radiological monitoring devices are quite variable in 
their accuracy of intensity measurements. Indeed, winds or redistributed fallout 
cause information inputs to reflect Incorrectly the status of Che radiological 
environment. Because of this uncontrollable reliability, decontamination de¬ 
cisions will have to be made despite some inaccuracy in information inputs. 




4. Comprehensibility 


Information inputs from external sensors must be collected (time, location, 
etc.) and presented in a comprehensible manner. Prestored data must be clear 
and written or presented in a manner readable by the user , e.g., decontamination 
planning guides should be readable to the local decision-maker in the £0C. 

D. Means for Obtaining, Storing, and Disseminating the Information 

Radiation information can be obtained from rr liation monitors. This informa¬ 
tion can be measured with instruments or visually observed. It may be forwarded to 
the EOC via telephone, radio, or hand-carried. Such information can be posted and, 
kept up to date on display boards, map books, or typed and hand-written reports. 

Field survey data can be handled in the same manner 

Prestored data are obtained by preattack u*ta acquisition and research. The 
kinds of prestored data obtained are determined by careful preattack assessment of 
postattack data needs. They can be stored as printed reports, maps, etc., or if 
retrie al machinery is available, they could be stored on punched cards, magnetic 
tapes, drums, discs, etc. 

E. Sumnary 

It is not the purpose in this appendix to prescribe precise']' 1 an information 
subsystem for decontamination. Whether or not a separate information subsystem for 
decontamination as opposed to, say, firefighting is needed is not answered in this 
ap;endix The answer to this question depends largely o. the organizational structure 
and system design of whatever higher e'helon systems are developed (e.g., Emergency 
Operating Systems). 

This appendix simply Indicates a basis for designing a decontamination informa¬ 
tion system. 
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V. CONCLUSIONS AND RECOMMENDATIONS 


A. Conclusions 

Some command and control structure will be necessary to conduct effective munici¬ 
pal decontamination . Whether or not a separate decontamination command and control 
system is necessary is not determined in this report--nor will it be determined with¬ 
out a far more extensive research effort. The decisions and activities related to 
decontamination do, however, require an elaborate information subsystem consisting 
of organized data files containing prestored (NFSS data, decontamination planning 
guides, etc.) data and postattack assessments (including system feedback). A good 
argument for a decontamination information subsystem could be presented on the basis 
of the complex network of information flow shown in Figure F-7. 

B. Recommendations 

It is recommended that decontamination be considered in the development and de¬ 
sign of any Postattack Emergency Operating System to be operational at the municipal 
or community level . Ic is further recommended that the systems requirements for con¬ 
ducting decontamination be analyzed more thoroughly and the necessary system functions 
and minimal organizational structure for conducting effective municipal decontamina¬ 
tion be more precisely identified. 
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